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AFOSR Grant No. 85 - 0130 


Starting Date: May 1, 1985; Date of this Report: May 30, 1988 

Abstract 

During this period I have extended the second-quantized derivation of 
quantum 1/f noise shown in the Second Annual Report to the general case of N 
particles present in the final state. I also have derived the quantum i/f cross 
correlations and the corresponding cross -correlation spectra, which are 
important for the calculation of quantum 1/f noise in Kinetic coefficients such as 
the mobility and the diffusion coefficient of the current carriers in solids. In 
order to better explain the foundations of quantum 1/f theory, I have given a 
derivation of the quantum i/f Schroedinger fields from quantum electrodynamics 
with the use of coherent states. Finally, I have given a direct derivation of the 
quantum 1/f effect in time and space. In terms of applications, a quantum 1/f 
noise study of MIS detectors was performed. Experimentally, with the 
collaboration of the group of Prof. A. van der Ziel, an excellent experimental 
verification of quantum 1/f theory was performed on semiconductor diodes, 
transistors and vacuum tubes, and a review article on the results of the 
experimental aplication and verification of my theory was published by A. van der 
Ziel in the Proceedings of IEEE in March 1988. v... 
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I. INTRODUCTION 

The present report will give a general description of the very extensive work 
performed in the reported period. This work was done along two major directions: 

1. Theory, which is presented in Part I of this Report. This had the goal to 
further develop the quantum 1/f theory and to clarify its foundations and 
starting points. This part has nine sections, an Appendix, and its own 
references. 

2. Application to MIS infrared detector structures, presented in Part II. This 
part is dedicated to a detailed quantum 1/f noise study of an important type of 
infrared detector, the metal - insulator - semiconductor structure. 

Part III lists the publications corresponding to this grant period. 

At the III Conference on Quantum 1/f Noise and 1/f Noise in Minneapolis, 
April 28-29, 1938, Prof. C.M. Van Vliet also presented a quantum 1/f noise 
derivation in the Van Hove limit. At present our attention is focussed on 
applications of the newly calculated cross - spectra of quantum 1/f noise. 
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PART I 


THEORY 







II. DERIVATION OF THE PAIR-CORRELATION FUNCTION FOR BOSONS 
Denoting by $(l?,t) and x(x.t) the single-particle wave functions of two 


scattered bosons, which take Into account both the Interaction with the 
scattering center of forces and the simultaneous Interaction with the photons, 
the quantum state of two bosons, both emerging from the same interaction process, 
can be written in the Heisenberg picture in the form 

|S> » (1/2) 1 / 2 jd 3 5jd 3 n^t,0) x (?r,0)/(5)l» + (^|0>. (2.1) 

where y(t) designates the Heisenberg field operator at a fixed time t = 0 and 
the single-particle wave functions have also been taken at t ■ 0. The two 
single-particle wave functions Included differ only through some external and 
Internal phase factors. Vectors are printed In boldface. The state In Eq. 

(2.1) Is written as a product of two single-particle wave functions, because we 
neglect all Interactions between the particles except for the quantum exchange 
symmetry between Identical particles. We have limited ourselves here to two 
outgoing particles for the sake of simplicity only; the general N-partlcle case 
Is treated In Sec. VIII. 

The operator of the equal-time pair-correlation 6 in space is 

0 « q» + (xi)q» + (x2)'p(x2)’*'(^i)- (2.2) 

This corresponds to a density autocorrelation function 6 . The presence of two- 
particle coordinates In the operator 0 does not mean that we are considering 
two-particle Interactions, It only means that the expectation value which we are 
calculating depends on the relative position of the particles. Using the well 









known commutation relations 


♦(x)* + (y) - ♦ + (y)*( x ) 3 «(x - y). 

(2.3) 

♦(x)*(y) - <|>(yMx) * 0, 

(2.3a) 

* + ( x )/(y) - ♦ + (y)/(x) 3 0, 

(2.3b) 


we obtain the matrix element: 

<S°|0|S°> 

3 (l/2)<0|Mn , )*(r)/(x 1 )/(x 2 )v(x 2 )^(x 1 )* + (|)^(n)|0> 

3 (1/2)[$U' - XiMn' - x 2 ) + <5(n' - xi)«(T' - S^)] 

[6(n - x 2 )«(5 - xj) + 6(5 - x 2 )6(iT - xj), (2.4) 

where |S°> is the state with well defined particle coordinates. 

The pair-correlation function is then 

A =« <S|0|S> » (1/2)< x *(x2)**(xiH(m)x(x 2) + X*(x2)#*(xi)^(x2>x(*l) 

+ X*(^i)4>*(x2)<>(xi)x(x2) + X*(xi)**Cr 2 )*(X2)x(xI)>. ( 2 * 5 ) 

1) If we assume that the wave functions <j>(x), X ( x> )» •••of different 
particles in the outgoing flux differ only through a general phase factor, we 
obtain from Eq. (2.5) 

<A> 3 2<U(xi) | 2 |$(x 2 ) | 2 >, 

which is similar to our previous result 2 2<|<^(t)<j>(t+-r) | 2 > 3 2< |<t>(ti)«|>(t 2 ) | 2 >, 
the only distinction being the use of the spatial coordinate along the beam 
instead of the time coordinate. As before 1 " 5 we write from Eq. (2.5) 

<A> = 2<U no (x 1 ) + Ei ^ i Br( x l)l 2 IW*2> + -i *V( X 2) l 2> 

3 2<(♦noI 2 + Ei U’Brl 2 ) 2 

+ 2<^* no ( x i)<t>no^ x 2) El < ^* 1 Br( x 2)^ 1 Br( x l) > + cc * 


( 2 . 6 ) 









Here we have separated the part with bremsstrahlung into the photon mode i, <^Br 
from the main part, $ no , without bremsstrahlung, and cc denotes the complex 
conjugate of the preceedlng term. Each <^ Br term has an Independent random 
phase; as we shall see In Sec. IV, this is the random In Itlal phase of each 
photon mode i, where 1 includes both the wave vector and the polarization. The 
averaging in Eqs. (2.6) and (2.6a) Is with respect to these phases, or, if the 
emitted photons are included Into the single-particle wave function, the average 
includes also the expectation value in the space of photons. Due to the uniform 
motion of the scattered particles, the outgoing radial coordinate and the time 
are equivalent; we will check this intuitive conjecture later in Sec. VI and we 
will prove by direct calculation of the correlation in space and in time that it 
is a very good approximation. 

We have denoted here by $ the stochastic Schroedinger field <j> used in our 
previous publications 1 " 5 . This field was not second-quantized, and not an 
operator as far as the charged particles were concerned. It was always called a 
field rather than a wave function, in order to emphasize the presence of random 
phase factors in its expression, which make it stochastic in nature, while the 
term wave function was reserved for pure states rather than mixture of states. 
The stochastic Schroedinger field description is equivalent with the density 
matrix description of quantum mixtures. The reason why still was an operator 
in the first paper 1 is that while the particles were not second-quantized in 
that paper, the electromagnetic field was. Therefore, in that paper n> was an 
operator only in the space of the photons. 

2) If, on the other hand, the single-particle wave functions $(x), x(x), ... 
also differ through mutually independent sets of random phases in their energy 
loss (<^Br» x^Br. »••) parts as we conclude in Sec. IV, the first and the last 
terms in Eq. (2.5) are constant and do not yield 1/f noise, but the two middle 
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terms give again the same result as In Eq. (2.6) 

<A> - <|*(x 1 ) | 2 ><|*(x 2 ) | 2 > + |<<J>*(x 1 )4>(x 2 )> I 2 

- (Unol 2+ S 1 U i Brl 2 ) 2 + l$*no( x l)*no( x 2) + W^l^Br^) l> I 2 

* (Unol 2 + s i I^Brl 2 ) 2 + Unol 4 + ( £ 1 l < > 1 BrI 2 ) 2 

+ ♦*no( x l)tno( x 2) Z 1 Br( x 2 )t 1 Br( X 1 ) + cc 

- <l<t*(xi) | 2 >< I<fr(x 2 ) | 2 > + <U(xi) | 2 U(x 2 ) | 2 >. (2.7) 

Both In the first and last form we notice that the first term Is constant, while 
the second is the familiar APSPSWF. Therefore, this Is also similar to the 
APSPSWF result obtained In our previous papers 1 " 5 , but this time the fractional 
spectrum will be two times smal1er f due (see Sec. VIII) to a 1/N factor with N*2; 
fermions have 1/(N-1). Here the average is with respect to the sets of random 
phases present in the energy loss (bremsstrahlung) contributions, or, if we 
Include the photons in second quantization into the final state (see Eq. 4.15), 
the average sign in Eq. (2.7) also includes the expectation value in the Hilbert 
space of the photons, which yields exactly the same result. The arguments have 
been omitted in some of the terms which do not actually depend on them according 
to Eq. (4.15) of Sec. IV.. Sec. VIII gives the APSPSWF result for any N. 

We conclude that both if the outgoing particle wave functions $(x), x(^)» 
..., differ only by an arbitrary general phase factor, and if they also have 
independent sets of phases in their bremsstrahlung energy loss parts, our 
previous method of calculating the autocorrelation function is justified on the 
basis of the present second-quantized calculation; see Sec. IV-V for details. 

In Sec. IV, we will determine the correct single-particle wave functions (or 
Schroedlnger fields) which lead to Eq. (2.7) and to the fractional spectrum 
2aA/fN for N outgoing particles (Sec. VIII) with and without the inclusion of 
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the photons Into the final state of the scattered particles. 

Why do we allow for so many possibilities ? Do we need to consider both 
cases with and without the photons included into the final state? Often In 
electrophysics one describes the quantum motion of an electron In an external 
time-dependent, electromagnetic field of force. In this most convenient 
description 7 energy Is not conserved, and the electromagnetic field energy Is 
not Included Into the Hamiltonian which is time-dependent. The state does not 
include the field either, and is a nonstationary electronic state. This 
corresponds to the treatment^ presented in this section, I.e., without the 
photons Included In the state (Sec. IV). An equivalent, more fundamental, 
approach Includes both the charged particles and the field energies with their 
Interaction Into the Hamiltonian. Then the photons must be present in the 
state, as we did before 1 . The two methods are equivalent, and give the same 
results, with no contradictions. We consider here both methods. 

Finally, why do we need to consider both the case of Identical single¬ 
particle wave functions (except for a general phase factor) and the case of wave 
functions which have different sets of phases in their energy-loss parts? 

Because both cases may occur in practice. Indeed, usually the incoming 
particles In a scattering experiment have a random shift in space or time which 
scrambles their energy-loss phases and eliminates cross terms as if these sets 
of random phases would be completely independent for different particles. 
However, in emission processes, such as a-decay, the single-particle wave 
functions are not shifted; they all start leaking out at t = 0 when the 
radioactive source was prepared by chemical separation. Therefore we must 
consider all cases in spite of the inconvenience. For instance, by writing 

explicitly ip(x) - 9noU) + £ fW*) a + i and X (x) = Xno (x) + £ x i Br( x )a + i, we 
1 i 

obtain from Eq. (2.5) in which the asterisks are interpreted as crosses 
Indicating hermltian conjugation, the same result as in Eq. (2.7), by 
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Interpreting the brackets as vacuum expectation values. 


III. DERIVATION OF THE PAIR-CORRELATION FUNCTION FOR FERMIONS 
In the case of fermions, the state of two scattered particles emerging from 
the same Interaction Is similar to Eq. (2.1) 

|S ss *> « (l/2)l/2 Jd3 € jd3^(5,0) x (n,0)/ s '(t)^ + s'(n)|0>, (3.1) 

but now the field operators satisfy anticommutation relations: 

*s(x)»V(y) + (?)*$(*) * 6 (* ■ y)«s,s'. (3.2) 

♦s(x)+s'(y) + 'I's'(y)+s( x > 3 °* (3.2a) 

♦ + s( *)♦%'(?) + ♦ + s'(y)» + s(x) 3 0* (3.2b) 

The operator of the palr-correlatlon function Is now written In the form 
0 - E S s' * + s (^ 1 )^ + s*(* 2 )♦s'( x 2 ^♦s( x l)• (3.3) 

Its expectation value In the state given by Eq. (3.1) is calculated with the 
help of the matrix elements 

<s hm\ 

- <0 l^n' )♦,(? ) ♦* (*1) (x 2) <&( x 2) ♦/ x 1) H» + (T)% + ("^) l 0> 

- C- 6 (n' -■*i)«U' - x 2 ) + «U* - x l)«(n* -~ x 2 )] 

C«Cn - X 2 ) 6 U - x i) - 5(7 - ^ 2 ) 6 (^ ’ x l)3; (3.4) 

- <01 ♦( n’) # f ( 7 ) **( X 1) +*( x 2) ^ x 2) * t ( X 1) ♦*( %) ♦*( n) 10> 

■ 6 (n* -72)6(5' - >4)5(7- x 2 )6(5 - xi); (3.5) 

<S°J0|S!> 

*r n if _ 

■ <01 ♦ (n') ^(?) ♦ J 1 (xi) ♦*(x 2 ) ^(x2) *^ xi) *+( 5 ) ^(n) 10> 

* 6 ( 5 * - X 2 ) 6 (n' - xi)6(5 - X 2 ) 5 (n - xj). (3.6) 

We also obtain three similar expectation values with all spins reversed. 

The spin-averaged pair-correlation function is then 









* (1/4)[|x(m)| 2 |*(x 2 )| 2 ■ X*(xi)#*( x 2)x(x2M x l) 

"X*( x 2)^*( x l)x( x l)^( x 2) + lx( x l) I 2 U( x 2 ) I 2 

+ lx(x 2 ) I 2 U(xl) | 2 + Jx(xi) I 2 |<t*(x 2 ) | 2 ] 

- (1/2)Cx*(x 2 )/(x 1 )^(x 1 ) x (x 2 ) + x*(xi)^*(x 2 )^(x 2 )x(x 1 )] 

- (1/4)[ x *(xi)/(x 2 )^(x 1 )x(x2) + x*(x2)+ # (xi)^(x 2 ) x (xi)]. (3.7) 

1) If we assume that the wave functions $(x) and x(x) differ only through a 
general phase factor, we obtain from Eq. (3.7) 

<A> = (1/2)< l<j>(x^) | 2 |*(xJ) | 2 > (3.7a) 

which corresponds to our earlier APSPSWF result 2 in the time domain, and which 
yields again the fractional spectrum 2aA/f found in our previous papers. (See 
Sec. IV for the definition of aA.) 

2) If, on the other hand, the single-particle wave functions also contain 
mutually independent sets of random phases In their energy loss parts, we obtain 
from Eq. (3.7) 

<A> =* <|*| 2 ><|*| 2 > - (1/2) |<$*(xi)$(xjT)> | 2 

3 (Unol 2 + z i l+ i Br I 2 ) 2 - (1/2) | U no I 4 + Ui U i Brl 2 ) 2 

+ t*no( x i)l>no(x 2 ) Si ^* i Br( x 2 )^ i Br(xi) + cc) 

* <|<M 2 > 2 - (1/2)<|MM)| 2 |^)| 2 >. (3.7b) 

Both in the first and the last form above we notice that the firstterm Is 
constant, while the second is the familiar APSPSWF. Therefore, this result is 
again similar to our previous APSPSWF result* -2 shown in Eq. (2.6) or (3.7a), 
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and yields the same fractional spectrum 2aA/f. However, at small distances we 
now obtain something similar to a Fermi hole due to the minus sign in Eq. 

(3.7b). The minus sign corresponds to a 180° phase shift or delay. This means 
that the fermions avoid each other at close range. Over large distances, l.e., 
at the low wave numbers significant for 1/f noise, the distribution and the 
spectral density are not affected. If we expllctly Include the photons into the 
final state and multiply <^Br b y a Photon creation operation a^j, we obtain 
again exactly the same result shown In Eq. (3.7b). This shows that although the 
phase shift caused by the Pauli principle at short distances was not included, 
our previous treatment 1 " 5 was correct because It started with exactly the same 
product of single-particle wave functions and therefore yielded the same final 
result as the present calculation. Whether this result should simply be 
interpreted In terms of frequency beats Is a philosophical question; I know It 
should. However, while the beats are caused by exactly the same Interference of 
each single-particle wave function with Itself which we had considered from the 
beginning, they occur because the single-particle wave functions are Included as 
symmetrized products into the many-particle wave function of the outgoing flux. 

The Fermi hole is absent In Eq. (3), because there only the antiparallel spins 
contribute. 

We finally conclude for fermions as well as for bosons, that both if the 

outgoing particle wave functions <fr(x), *oT).differ only by an arbitrary 

general phase factor and If they also have independent sets of phases in their 
bremsstrahlung energy loss parts, our previous method of calculating the 
autocorrelation function Is justified on the basis of the present second- 
quantized calculation; see Sec. IV-V for details. All claims to the contrary 
which have been recently voiced, are therefore wrong. We now need to insert the 
correct single-particle wave functions explicitly into Eqs. (2.5) and (3.7). 
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Before we do this direct calculation, we briefly discuss the connection 
between density and current density autocorrelations. The current density 
operator In second quantization Is 

J « Cli/21mh + (5t)(M)*(x). (5.6) 

where now *(x) and \|» + (x) are particle field operators, while the same notation 
was used earlier 1 to designate single particle wave functions which are 
operators In the Hilbert space of photons, because the current-carrying 
particles were not second-quantized, while the emitted photons present In the 
final state were quantized. The operator of the current correlation which 
replaces Eq. (2.2) Is 

Oj - (V2m)V(xi)* + (? 2 )(^2)*(^2)(Ml)*( J l)» (5.7) 

while In Eqs. (2.5) and (3.7) we have to make the substitutions 

♦ (*l)'- CR/21m)($!-$i)$(?i) - (HK/mUtfi) 

♦(*2) + tt/21m)($ 2 -$2H(*2) » (HK/mM* 2 ) (5.8) 

and similar for xU)» while and x* are left formally unaffected. Here we 
have used the form of $(x) and x(x) given by Eq. (5.1) or (5.2) which allow only 
for very small momentum deviations q « K caused by infraquanta, which can be 
always neglected. Therefore, Eqs. (2.5) and (3.7) will only acquire a general 
factor (HK/m) 2 . We conclude that in general 

<S|0j|S> = (T>K/m) 2 <S |0 |S>, (5.9) 

l.e., that the current density autocorrelation Is proportional to the density 
autocorrelation If the single-particle states are close to momentum eigenstates. 
The corresponding spectral densities are also proportional, with the same 
proportionality factor (HK/m) 2 . The fractional autocorrelations, however are 
the same for the current, as for the density, because they have the squared 












average value In the denominator. The same equality holds for the corresponding 
spectra. 

Finally, for single-particle states which are superpositions of states close 
to momentum eigenstates, the above proportionality remains approximately valid. 
If we replace the proportionality factor with <(fiK/m)2>, where the average Is 
over the distribution of K . 








IV. DERIVATION OF THE SINGLE-PARTICLE WAVE FUNCTIONS 


The derivation of the Schroedinger field will be performed with the 
help of the Green's function method similar to the method used in a earlier 
calculation by Kroll and Watson 7 , extended to the case of interaction with 
all electromagnetic modes of the universe*. 

It is most convenient to describe the electromagnetic field in terms 
of plane waves. The vector potential is taken in the radiation gauge as 

A(r,t) = 2 »,i(h l d l /L J 6)k) 1/2 Uk,itak, i(t)e lk< r + a**, i(t)e- 1 k * r ]. (4.1) 

The polarization vectors u»,i and Uk,a are mutually orthogonal unit vectors 
perpendicular to k. 

The Schroedinger equation for an electron moving in a vector poten¬ 
tial A and scattering potential V is 

(l/2m)[-iliV - eA/c]»<4 + V0 = m. (4.2) 

A dot has been used to indicate the time derivative. The electromagnetic 
field is treated as a classical field at this point. In order to eliminate 
the A* term from Eq. (4.2), we write 

V - expn-i/hJJMeVZmc’jA’dt’]# . 

Thus, Eq. (4.2) is reduced to 

((-■nV2m)V* + (ieft/mc)A*V + V]<*> = ifi<j> . (4.3) 

It is convenient to consider first the influence of a single electromagnetic 
mode, i.e. a single term from Eq. (4.1). Therefore, we take A » acos( oit+r) , 
where y is an initial phase constant, and we treat V<X> as a perturbation 
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source term. The solution for Eq. (4.2) is an incoming plane wave plus 
scattered waves, given by the integral equation 

#kjr,t) = 4>* c - XdJx’p dt’GV(r’)<l> k c (r\t'). (4.4) 

Here is the solution of the homogeneous equation, i.e. with V = 0, and 
can be written in the form 

= e‘»*rexp[-(ifi/2m)J‘(k* 0 - 2ek 0 -A/ch)dt]. (4.5) 

G is the Green's function which satisfies the equation 

[(-h*/2m)V * + (iehA-V/me) - ihd/9 t]G = $(r-r')£(t-t'). (4.6) 

Given A = acos(<ut + r)» G can be found to be 

G = [i/(2x) J R]Jd 1 ke 1 k (r-r’ )exp{-(ifi/2m) [k 2 t - 2ek-asin(<ut + r)/hca>]} 
x exp[(ift/2m)(k*t' - 2ek-asin(<ut( + r)/h ca], (4.7) 

In the first Bom approximation we set <J> k fl (r 1 ,t') = <£» o (r* ,t’) in the in¬ 
tegral present in Eq. (4.4) and obtain for the scattered wave 

= [i/(2*)»!i]Jd*x , J , » dt’V(r’)Jd 3 ke•*<*-*' >exp{-(ih/2m) 

[k J t - 2ek*asin(<ut + r)/hc<u]) x exp[(ifi/2m)(h* - k*o)t'] 
x {exp(ie(k# - k)asin(6Jt' + r)/2mc<u]e* k e ' r ’) • 

Using the relation 

(4.8) 

where J«(/9) is the n th order Bessel function, we expand the expression con¬ 
tained in curly brackets in Fourier series. Then Eq. (4.4) takes the form 
<Dk # (r,t) - ** 8 (r,t) = (-i/(2*) Wd*x’.pdt'V(r')Jd*ke>*<r-' , > 
x exp{-(ifi/2ra)[k l t - 2ekasin(<yt+r)/ftco>]} 

x exp((ifi/2m)(k* - k*)t’][Z (4.9) 

After performing the integration over t’ we use a contour integration method 
for k. Then Eq. (4.9) is reduced to 

$* e <r,t) - <f>k c (r,t) = [-m/(2jr)fi 1 ]Z J«(;3)e‘ »*[e>*(>>'/r] 
x exp{ift[k*(n)t-2ek(n)asin(cjt +r)/hc<w]/2m} x Jd , x’e 1 r 'V(r')e l k « r ', 
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$ - -e(fik* - tiko) -a/mctia) = -eQ-a/mchai, 


(4.11) 


and Q is the momentum transfer. In Eq. (4.10) k(n) is defined by 

(hk(n))*/2m = Cfck«,)*/2m - nho>; k, = k(n) = k(n)r/r. (4.12) 

The total scattered wave can be written as 

= [-ra/(27r)f»*r]Z «.- J »e ,k < > > r exp{-iTi[k J (n)t -2ek(n) •asin(<»>t+r)/hc6)]/2m) 
x V k(l)>k J,(^)e‘«^ (4.13) 

where Vk(a),k = Je* lt(,)r, V(r , )e* k o ,r ‘d , x’ is the scattering matrix ele¬ 
ment calculated without consideration of the interaction with the 
electromagnetic field oscillators, bnenltutici t 0 Jl foldoyullthn*pkiu f &r exe X i 
So far the electromagnetic field has not been quantized and was con¬ 
sidered as a classical field. We are interested in the corresponding ex¬ 
pression of the scattered single-particle wave function when the electromag¬ 
netic oscillators are quantized. Therefore we first linearize Eq. (4.13) 
with respect to the electromagnetic potential wherever a dependence on § is 
present: 

= [-m/(2x)tt I rJVk,k e 1 «*- 1 '‘'' , exp[ieK-asin(<Jt + r)/cma>] 

{1 + e i(uiMr*/)|3/2 - e - 1 >>9/2} 

= (-ra/(2»)fi , r]Vi,k e‘* r - ,l ‘/*exp[l + eiK-asin(a»t + r)/cm6>] 

{1 + e i <«/»-, fneQ*a/2mdfi<i) - e-« <‘"»-<r*/) e Q. a /2mrfi6>}. (4.14) 

Since k = (v/c)q<<q, the r-depedence of the elfctromagnetic potential can be 
neglected, as it is, in exp [ieK-asin(cut + r)/cm<u]. 

Here we have introduced the notations E = (fiK)*/2m and K = k(0). In the 

last form we will neglect the terra with sin (cut + r) tfcat corresponds 
part ef-.’tk* 

to a\coherent quantum 1/f noise contribution which has been considered else¬ 
where before"* 11 , because this term can be considered constant of negli¬ 
gible magnitude in the calculation of equal time spatial correlations. We 
also neglected the small difference between k(l) = K - ojK/vK and K in 0 and 
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v k(1),k 0 • We conclude that the quantization of the electromagnetic field 
transforms <|» s into an operator in the photon Hilbert space 

4>(r,t) = (C/r)e'“ r *> *»/*{1 - 2 k , ib*(k,l)expi(-«jt + qr - r i)a k ,i 

+ 2 k, ib(k,l)expi(<ut - qr + 7 i)a f k,i). (4.15) 

This Is the form of the single-particle wave functions which are used In the 
present paper with or without the photons Included. Here we have Introduced 
b(k ,1) ■ (1/2)0 and the constant C which designates the factor In front, and a 
sum which Includes all electromagnetic modes with annihilation operators a£ 

In the space of the electron states $(r,t) is just a single-particle wave 
function. We have denoted by q the small decrease In the particle momentum 
required by Eq. (4.12). We have q * (K/E)*Hu » ck/v * w/v, with 

(L/2*) J -4>r-Z i<| b(k,l) | *>k*dk = e*Q*aV4m«c*h*<w*-2k*d<i>/3c = aAdm/m, (4.16) 


where a - e*/4«fic = 1/137, A = (2Q*/3mn J c*), and where < > is an angular 
average. We have considered the spontaneous emission caused by vacuum fluc¬ 
tuations only, yielding for the amplitude Cy, 




The annihilation part is included in Eq. (4.15), but does not contribute to 
the quantum 1/f noise on the background of the electromagnetic vacuum. If 
the calculation is performed on the thermal radiation background, however, 
we get a white noise contribution added to the quantum 1/f noise which 
remains the same*. Here we have performed the transition from just one 
electromagnetic mode to the general case with all electromagnetic modes ad- 
hoc, but in our previous paper 8 this transition was presented in detail. 
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In Eq. (4.15) we notice the presence of the random phases n which 
where introduced as initial phases of the electromagnetic oscillators, and 
which are independent for each electromagnetic mode i (of given polarization 
1 and wave vector k) of the universe. Since the various scattered particles 
are independent of each other and of the electromagnetic modes, we have to 
consider the set of random phases 7i different and independent in the wave 
function of each particle, als« Me disco-ssCc*. before E^.(J./p)aW after fy.(AUf). 

V. CALCULATION OF THE PAIR-CORRELATION FUNCTIONS AND 
OF THE CURRENT AUTOCORRELATION 

In Sec. IV we have determined the single-particle wave functions 
in the form, valid at t = 0, 

tf(x,t) = (C/x)exp(iKx)[1 + Z *,ib(k,l)exp(-iqx)], (5.1) 

or, if the photons emitted in the bremsstrahlung process are included in 
second quantization, 

tf(x,t) = (C/x)exp(iKx)[l + Z k ,ib(kl)exp(-iqx)a* k , ij, (5.2) 

Substituting these expressions into the calculated expectation values in Eq. 
(3.7), we obtain for the case of fermions 

A(xj, x 2 ) = (C*C/xiXa)*{1/2 + £ k .i|b(k,l)| *[2-cosq(x x - x»)] 

+Z kk'ii'l b(k,l) | *| b(kM') | *[1 -(l/2)cos(q -q’)(xi - x 2 )]}. (5.3) 

Note that if Eq. (5.2) is substituted, the expectation on the electromag- 







netic vacuum state is calculated. If the thermal radiation background is 
also taken into account, a small white noise term is obtained in addition to 
the 1/f noise*. The constant part is the squared expectation value of the 
concentration of particles. Dividing the variable part by this constant 
part,we obtain the fractional spectral density of the particle concentration 
n, or current density j (with L J being the volume of the normalization box 
considered) 

S.(k)dk/n* = Sj(k)dk/j* = 22 i«| b(k,l) | »>>(L/2*) 3 4xk*dk 
/{I + 42 k.i | b(k,l) | * + 22 kk* i i>| b(k,l) | *| b(k’,l’) | *} 

= 22 i«| b(k,l)J *>>(k , /2w 1 )L 3 dk = 2oAdk/k = Sj(f)df/j* = 2aAdf/f, (5.4) 

which is in agreement with our previous results, and which also includes a 
180* phase shift due to the exclusion principle which is important only at 
short distances between the particles. In the final form we have trans¬ 
formed to the frequency f. 

The expression of b(k,l) used in Eq. (5.4) was derived in Sec. IV. 

It is the well known expression of the bremsstrahlung amplitude in any scat¬ 
tering process 

b(k,l) = - ( e/riKa) (h/ai. J ) 1 /*A pAi > (5.4a) 

where Ai is the polarization vector of the mode, and Ap is the momentum 
change of the particles in the scattering process considered. In Eq. (5.4) 
an average << >> over the angular variables in k and a summation over 

photon polarizations were performed, leading to the usual definition of the 
quantum 1/f noise coefficient, also known as the nonrelativistic form of the 
infrared exponent defined in quantum electrodynamics 

oA = (2a/3jr)(Ap/mc)*, (5.4b) 

with a being Somnerfeld's fine structure constant a = e*/ftc = 1/137. 
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In the case of bosons we substitute Eq. (5.1) or (5.2) into Eq. 


(2.5b) and obtain 


A(xi.xj) = 2 + 22 k, i | b(k,l) | *[1 ♦ cosq(xi - x»)] 

+ 2 *,i2 k*. i • | b(k,l) | *| b(kM') | *[1 + cos(q - q»)(xi - x*)]; 

S a (k)dk/(n) 1 = Sj(k)dk/(j ) 2 = 2 i«| b(k,l) I ’»<L/2jr)* 4xk*dk 
= aAdkA = aAdf/f. 


Here at short distances we notice an increase of A(xi»Xi). The fractional 
spectral density is reduced by factor of two compared to the case of fer¬ 
mions. Both results generalize our previous results to the case of short 
distances and prove the existence of the conventional quantum 1/f effect in 
second quantization, even if the photons are included into the final state 
(which is to be avoided if the electromagnetic field is included in the 
Hamiltonian as an external field, as we did in earlier papers). 

The transition from the wave number spectrum to a frequency spectrum 
in the last form Eqs. (5.4) and (5.5) is based on the equivalence of the 
density distribution along the outgoing rays with the time-dependence in a 
fixed point along the flux of particles which all move with the speed v. 

The coordinate x along the ray and the time t are in this case equivalent 
variables, and therefore the wave number k and the frequency f are also 
equivalent variables for the spectral density. We can therefore write 2xf = 
vq = ck and S(k)dk = S(f)df, which justifies the last form of Eqs. (5.4) and 
(5.5). In order to check the validity of this procedure, we will now per¬ 
form a direct calculation of the pair-correlation function in time and 










VI. DERIVATION OF THE PAIR-CORRELATION FUNCTION IN TIME AND SPACE 


Working in the Heisenberg representation as above, we can generalize 
the derivation presented above by including two different times in the 
operator of the pair-correlation, although this requires more calculation. 

To simplify the integrals, we consider again the state of two outgoing fer¬ 
mions very far from the place where they have independently suffered the 
same interaction, so that the outgoing spherical waves can be approximated 
by plane waves. Starting therefore with plane waves similar to the spheri¬ 
cal waves used in Eq. (3.1), with q = qK/K = ckK/Kv and q’ = ck’K/Kv, we ob¬ 
tain again the pair-correlation function in the form 

A(xi, tr, xa, ta) = l/2.fd*Sd’7jdi£'d37?» 

{exp(-i& 7 ») + 2 ,s' )exp[-i(K-q’ ]aa>,>) 

(exp(-iK£») + 2 »..b*(k,s)exp[-i(K-q)«']a k ,} 

(exp(iK£) + 2 *,.b<k,s)exp[i(K-q)£]a**.} 

(exp(iK^) + 2 k*,*-^(k’,s , )exp[i(K-q , )7]a**»,>) 

(1/4)2 ss*<S°s3• | O f| + <y 0^1 I( (t , ; . » 2 | S*ss- >. (6.1) 

The operator of the pair-correlation function now contains two consecutive 
times. The creation and annihilation operators for particles obey anticom- 
mutation relations similar to Eqs. (3.2)-(3.2b). Using these operators, the 
field operators can be expanded in terms of plane waves 

0s(r,t) = V-w*2 p expi(pr-cut)c»s, 0 4 s(r,t) = V-»/*2 p expi(-pr+<ut)c**s(6.2) 

which also contain their time dependence, as needed for the operator of the 
pair correlation function in the Heisenberg representation. This operator 
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has an expectation value which can be written f for spin up only, in the form 

<S# « I C V Xl * tl ’ x *» tl) I S V 

= < 0 | )^*(xi,ti)^(xa,tj)^(xi,tj)0(xi,ti)0*(5)0*(7) I 0 > 

= (1/V«)Z ■ expi[«* 7 ' +n’£' -/*’xi +<u(/i’)ti -*'xa +o>(i/’)ta 

+/oca -<u(^i)ta +ioci -<a<i/)ti -m >7 -n£] 

< 0 | c n <Cii'C y *, c£, c v c*«c*»| 0 > 

= (M/hjHt-Ha-^expfOM/ZfiJKxt-^) */ti +(x,- 7 )*/ta 
-(xi- 7 ')*/ti -(xi-fMVtsl) 

+ (M/h)•1 1 ■ J t 1 ■ J exp{(iM/2fi)[(X 1 - 7 ) */ti + (xa-£) */ta 
-(xi-f')*/ti -(xa- 7 ')Vta]) 

-(M/h)*ti->ta- 3 exp{(iM/2fi)[(xi-5) , /ti +(xa- 7 )*/ta 
-(xi-f’)Vti 

-(M/h)«ti- 3 ta- , exp{(iM/ 2 fi)[{xi- 7 ) */ti +(xa-f) */ta 

-(xi- 7 ')*/ti -(x,-f')Vt,]). (6.3) 

Substituting into Eq. (6.1), and integrating with respect to £, 7 , %' and 
if , we obtain for the part with spin up only 

8 Ajj(xi, ti,*xa, ta) 

= <0(1 +Z k'/SMk’ )expi[q'xi -iiq'tK -q’/2)ti/M]ak*> 

(1 +Z kb*(k)expi[qxa -Hq(K -q/2)ta/M]ak) 

(1 +Z kb(k)expi[qxi +hq(K -q^tj/MJa*a) 

(1 +Z k*>9(k*)expi[q’xa +1iq'.(K -q’/2)t,/M]a*k«}0> 

= 1 +Z k-| b(k’)| J expi(qt(xi -x,) -ftq’.(K -q’/2)(ti -t,)/M] 

+Z k| b(k) | *expi[q(xa -x.) --hq(K -q/2)(t, -t,)/M] 

+ Z k,k- I b(k)b(k’)| , expi[(q -q’)(xa -Xi) -fiK(q-q’)(ta -ti)/M 
♦T»(q* -q* 1 ) (t a -ti)/2M]. 




(6.4) 









The first form of tuts result proves the presence, with t f #r y of 
four single-particle functions (APSP5WF combination) in A^. 

In a similar way the expectation values corresponding to Eqs. (3.5) and 
(3.6) are calculated and integrated with respect to £, n, €’ and 77* . Put¬ 
ting together the terms for all spin orientations as before, we finally ob¬ 
tain the desired pair correlation function in the form 

A(xi ,ti ;x*,t*) = 1/2 + 2fc,i| b(k,l) | M2 - cosq[xi-xa -v'(ti-ta)]} 
+2 fcfc'i 1 • | b(k,l) | 1 1 b(k* ,1) I Ml -<l/2)cos(q-q*)-txi-Xa -v* ’(ti-ta)]);(6.5) 
v’= v( 1 -q/2K) = v(l-e/4E) = v; v'* = v[l -(q+q’ )/2K] = v[l ~(e+e’)/4E] = v. 

The approximations v* = v and v’ ’ = v are justified because the soft photon 
energy e = 4-10" ls eV for 1Hz is negligible compared to the energy E of the 
particles which may be of the order of leV. This result directly proves our 
heuristic generalization of Eq. (5.4) to the case of different times in the 
pair-correlation function. The extension to the case of bosons is trivial. 

VII. CROSS-CORRELATIONS AND SPECTRA 

In practical calculations of quantum 1/f noise in condensed matter 
and electronic devices we often need the cross-correlation of the differen¬ 
tial scattering cross sections at different angles. This is in fact the 
cross-correlation of the outgoing current densities scattered into different 
directions, with different incoming wave vectors Ki and Ka considered. If 
we restrict ourselves again to spatial correlations for simplicity, we need 
to examine the cross-correlation of the probability densities in the scat¬ 
tered wave at various distances from the scattering center and in various 
directions. 

For bosons we start from Eq. (2.5) and consider ^>(x) and x(x) as in 
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Eq. (5.2), but different from other; 

<£(x,t) = (C/x)exp(iKix)[1 + Z n,ib(k)exp(-iqix)a*k,i], (7.1) 

*(x,t) = (C/x)exp(iKax)[1 + Z k,i/9(k)exp(-iqax)a*k,i], (7.2) 

where now b corresponds to the momentun change K’ - Ki, and 0 corresponds to 
the momentun change K" - Kj of the second particle, in addition of having a 
independent random phase different from that in b, for each k, due to the 
different initial phases of the electromagnetic oscillators registered by 
the second particle. As in Eq. (2.5), we obtain now 

A(Ki,Ka,K’,K";xi,xa) = <S| 0| S> = (1/2) (| *(x,) | *| tf(x.) | * 

♦ X*(xa)4*(x,)x(xa)<t>(xi)+xHx,)<t*(xi)z(xi)<f>(xi) +1 2(xi) | *| <£(xi) | *} 

= [1 +Z k.1 I d(k,l) I *][1 + Z K.I I b(k,l) I 1] 

♦ (l/2)exp[iK*(xa-xi) ♦ iK“(xi-x*)][l + Z k,i| fl(k,l) | 2 exp-iq(xa-Xi)) 

11 ♦ Z k, i I b(k,l) | 2 expiq(xa-xj) ] + 11**—► 2] 

= [1 + Z k.i | d(k,l) | 2 ][1 + Z k.i | b(k,l) | 2 ] + cos(K' -K“)(xi-xa) 

+ Zk,i| d(k,l)| *cos[(K' -K") (xj-xi) + q(xi-xa) ] 

♦ Z k, i | b(k,l) | 2 cos[ (K’ -K")(xa-xi) - q(xi-xa)] 

+ Z kk-n- | d(k,l) | Mb(k’,l')l *cos(K’-K u +q , -q)(x 1 -xa). (7.3) 

Taking into account that the momentum changes caused by infraquanta are very 
small (q<<Ki-Ka)» we now smooth this rapidly oscillating correlation by 
averaging over short distances of the order l/(Ki-Ka) and get zero, except 
for Ki = Ka» which still allows each of the incoming directions and outgoing 
directions to be different 

A(Ki,Ka,K',K";x l ,xa) = <S | 01 S> S [ 1+Z k, i | /9(k,l) | *][1+Z k.i | b(k,l) | 2 ] 


+ (1 + Z k, it | /9(k,l)) | * ♦ | b(k,l) | 2 ]cos(q(xi-xa 
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♦ 2 kk'.ii’l y8(k,l)| *|b(kM')| *cos(q-q')(xi-x,)^ ,«* (7.4) 

This result indicates that only particles of the same energy have quantum 
1/f noise cross-correlations, and groups of different energy yield indepen¬ 
dent quantua 1/f noise, as was anticipated empirically by Kleinpenning 10 . 

Our result in Eq. (7.4) differs somewhat from a heuristic generalization of 
the basic quantum 1/f noise formula used in a previous paper®. 

For fermions we start from Eq. (3.7) and obtain 

A(Ki,Kj,K’ ,K";xi,Xi) = (1/2) [ | *(xi) | *| *<x«) | * + (l-->2)] 
-(l/4)[**(xi)tf*(xi)2(x,)4(xi) + (l«--2)] 

= [1 +Z k.i | d(k,l) | *][1 + Z k .i | b(k,l) | *] - (l/2)cos[(K'-K*)(xi-xi)] 
-(1/2)Z k.id /9(k,l)| * + |b(k,l)| *]cos[(K'-K")(x»-xi) +q(xi-x*)] 

-<1/2)Z kk*i■• I d(k,l) | Mb(k>,l’)| *cos(K‘-K"+q’-q)(xi-x*). (7.5) 

Applying again the smoothing operation we obtain now 

A(K l ,K 2 ,K’,K*’;xi,X 2 ) = [1+Z k.l | /9(k,l) | *][1+Z k.l | b(k,l) | *] 

-a/2){l + Z k,i[ | /9(k,l)) | » + | b(k,l) | *]cos[q(xi-xa)]}$r 

-(1/2)Z kk*, 11 * | /9(k,l) | *|b(kM»)| acos(q-q’)(x»-X2)5r (7.6) 

For elastic scattering we have K‘ = Kj and K" and Kg, provided the scattering 
center is fixed or very massive (e.g. a crystal). We Include the infrared 
radiative corrections this time. 

Using Eq. (4.16), the angular integrations can be performed in Eq. (7.4). 
We also apply the relations dk/k * df/f and k/k 0 - f/f 0 and obtain 
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AbtKi.Kg.K'. K",x lf x 2 ) - Cl + aAj j F (f/f 0 )° A ldf/f][1 + aA 2 J W 


+ 11 + aj F [A 1 (f/f 0 ) aA l + A 2 (f/fo) aA2 ]cosq(xi-x 2 )(df/f) 


(7.7) 


+ a^Aj, j F (f/f 0 ) aA l(df/f) /(fVfoJ^^Vf'Jcostq-q'Jtxi-XgJjfi^^,,. 

^0 *0 

with Ai * 2fl 2 (K , -Ki) 2 /3«m 2 c 2 and A 2 * 2fT 2 (K"-K 2 ) 2 /3nm 2 c 2 for bosons. The lower 
frequency limit f 0 Is the resolution limit given by the reciprocal duration of 
the noise measurement. The upper limit F Is very high, close to the energy of. 
the scattered particles, divided by Planck's constant. From Eq. (7.6) we obtain 

F F 

A f (Ki,K 2 ,K', K-.Xi.Xg) - Cl + oA x J (f/fo^ldf/fHl + a A 2 J (f/fo^df/f] 


-(1/2)11 + a/ CA^f/fo)®^ + A 2 (f/f 0 ) aA2 ]cosq(x 1 -x 2 )(df/f)j5 K 


(7.8) 


2 /2)A^A 2 J F (f/f 0 ) aA l(df/f) J F (f , /f 0 )° A2 (df , /f , )C0S(q-q , )(x 1 -x 2 ))« Kl>lcl 
fa *o 


for fermions. The Infrared radiative correction factors (f/fo)^ were 

calculated elsewhere®and can be set » 1 for all practical cases. Neglecting the 

"noise of noise" term proportional to a 2 , and using the Wiener-Khintchlne 

theorem, we obtain the corresponding cross-spectral densities of the scattering 

+ ♦ ♦ ♦ 
rate fluctuations aw* ♦ > leading from Ki to K , and aw* ♦ > from K 2 to K 
Kj,K k 2 .k 

S M (K lt K 2 XX') b - (a/2f)[A 1 (f/f 0 ) aA UA 2 (f/f 0 ) oA2 3<w*^* l ><w* 2t * ll >5 K , tK) , (7.9) 

for bosons, the.factor In front being a/N^f in. general (N»2 here), and 

W^l.^.K',K") f - (a/f)[A 1 (f/f 0 ) aA l + A 2 (f/f 0 ) aA 2]< w+ >< WJ >6 (7.10) 

Kj,K K 2 ,K * »* 

for fermions, the factor in front being a/(N K »-l)f now (see Sec. VIII). 
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VIII. PAIR-CORRELATION FUNCTION FOR N OUTGOING PARTICLES 


So far we have considered only the simplest case of two 
particles present in the outgoing state, which is the minimal 
number that allows for correlations between different particles 
to be defined. Quantum 1/f noise is the manifestation of long- 
range correlations between particles, causing them to be bunched 
and to be subject to superpoissonian statistics. Therefore, we 
are interested in these correlations which reveal themselves in 
the pair-correlation function, also known as the two-particle 
distribution function. 

In general, there will be more than two particles in the 
outgoing state, which are observed in order to determine the 
scattered current. In calculating the pair-correlation function, 
we should therefore consider the general case in which N scat¬ 
tered particles are simultaneously under observation in the out¬ 
going state during the measurement of the cross section. This 
is what we will try to do in tis section. 

Similar to the case of two particles in Eq. (2.1). the 
state of N scattered particles can be written in the form 


| S> = (N!)-»/*n | 0> = II ijd^i 4i(£i) | S«>. 


( 8 . 1 ) 












Here we have denoted by 4>(£i) the exact single-particle wave function of 


the particle i, where both the interaction with the scattering center and 

mriik I 

the interaction'the electromagnetic field inodes are taken into account. For 
bosons the operator of the pair-correlation is given by Eq. (3.3), and by 
using the permutation relations Eqs. (3.2)-(3.2b), we obtain 


N! <S° | 0| S*> = <0) M ?*)...M?i)tf + (Xi)tf*(xj)iMx 2 )tf(xi)tf*(fi )...*»♦(£») | 0> 

= <0 -X1 )vA -x 2 ^^ 1 ) 5(?,-X 1 ) 5(? -"X 1 ) I 0> 

= Z,' z' div*-X i)8(vm 2 )<S(£ b -xi)$(£.-x 2 )Z n' n |) 

¥ T ' - ' v d 

= Z Z *(>?A-Xi)$(»A-x 2 )*(£ b -Xi)$(£.-x 2 )Z n d{T?i-e .), (8.2) 

A A 


where the sun over permutations runs over all permutations of the N-2 in¬ 
dices j. This allows us to calculate the complete matrix element 


<S| 0| S> = [1/N(N-1))£, jB |i ,Jd»^Jd*^Xd»fiard*«i 

^*(^)^*(^)^ ,( e , i ) ^.(^. , <y<^-xi)<y ( ^-x 2 , ^ ( £.-xi ) «y(£.-x J ) 

= [l/N(N-l)]Z^j^A *(xj)^«(xi)^i(xx)d.(x*) 

= [1 ♦ Z k.i | b(k,l) | *]"*{N(N-1) + [(N-1)/N][2Z -iZ k.i | b.(k,l) | * 
N 

+ 2Z ..i Z k.i| b.(k,l) | *cosq(xi-xi)]} 

= [1 + Z k.i| b(k,l) | *]■*{N(N-1) + 2(N-1)[Z k.i I b(k,l) | * 

♦ Z k.i | b(k,1) | *cosq{xi-Xj)]). (8.3) 


Dividing the term dependent on xi-xj by the constant term, and neglecting 
terms of higher order than | b(k,l) | *, we obtain the fractional spectrum 
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density of the quantum 1/f fluctuations in the concentration c of particles, 
in the current density j, and in the physical cross section a 


Sj(f)/<j>* = 2aA/Nf. 


In Eq. (8.3) we have used the form of the single-particle wave functions 
given by Eq. (5.1) or (5.2), and in Eq. (8.4) we have used the expression of 
the bremsstrahlung amplitude given by Eq. (4.16). Our calculation in Eq. 
(8.3) and the result in Eq. (8.4) show how the independent sets of phases in 
the bremsstrahlung energy loss part of the single-particle wave functions 
lead to independent additive quantum 1/f spectral contributions from each of 
the scattered particles, and to the 1/N factor in the fractional spectral 
density shown in Eq. (8.4). This result and the foregoing calculation are 
in agreement with our previous results 1 * 5 and calculations, and perfectly 
justify our previous use of single-particle wave functions with subsequent 
introduction of the 1/N factor on the basis of the independence of 1/f noise 
contributions from different carriers. 

In the case of fermions the calculation is similar, except for the 
use of anticomnutators for the fermion field operators. In order to em¬ 
phasize the independence of our results on the representation used, and to 
show directly how the calculations presented in this paper can be performed 
without second quantization, we give here the direct calculation in terms of 
a Slater determinant for the state of N scattered fermions 


^u(ri) <f> u(rj).tfu(rit) 

^ia(ri) tfii(rj). <f> u(rn) 


*(ri...r>) = (N!)* 1 





m 

m 

n 

i 


I 

?’*/< 

m 

a 

• 

a» 

'!;Jv 

1 


i 


^m(ri) <bin(rj).^iN(rn) I (8.5) 
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Here r combines the position vector x and the spin variable s for each of 
the particles. The pair-correlation function is obtained by integrating 
with respect to the coordinates of all but two of the fermions 

A(ri,r 2 > = Jd J rj.. .d 3 r»< | ... m(ri...m) | J >. (8.6) 

Here the integrals also include summations over the spins, and the expecta¬ 
tion value is with respect to the phases present in the brerasstrahlung parts 
of the wave functions and, if the emitted photons are included in the final 
state, the expectation value is also done over the vacuum of the photons. 
Assisting orthonormality of the functions 4> we obtain 

H 

A(n.ra) = [1/N(N-1)]Z | tf.(n)*.(r,) - tf.(r*)*,(n) | *>. (8.7) 

To display the spin variables explicitly, we write - Zm(x) | s> and get 

Hi 

A,,.(xi,xi) = [1/N(N-1)]Z m,».i<[Xm*(Xi)<S\ Z,*(xt)<a' I (8.8) 

- *.*(X 2 )<S’ | *o«(Xl)<s| ] [*.(Xl)|S> 2.(Xa)|s'> - X.(X 2 )|S'> Z.(xi)fs> ]> 

Hi 

= [1/N(N-1)]Z I Xm(xi) I 1 I X.(x a ) I * + I x.(x 2 ) I *| *,(X!) | * 

- *.«(Xi)*.(Xa) | <s| 3* > 1 *Xi t {Xi)Xn{Xi) - (Xi«-•* X2 ) ] >. 

Here the simbol (xi+--»x 2 ) designates the immediately proceeding term. Con¬ 
sidering all spin orientations, we obtain 

H/i 

A(xi,x>) = V \ = 11/ N(N-l)]Z.,,.i<(4| Z.| »| z.| > 








■ Xm*( x l)x*n( x 2)xm( x 2)xn( x l) " < X 1 * * *2)3* ( form 1 ) 

« [|C| 4 /x 2 N(N-l)3lN 2 Cl+Ejr t i|b(k,1)| 2 ]2 
N/2 - 0 

-2 £ exp[1K m (xi-x 2 )]Cl + ijTj |b m (k,l) | 2 exp[iq m (xx-X 2 )] 
m-1 
N/2 

•E expC-1K n (xi-x 2 )]Cl + Eg* j • l b n( k * .1') I 2 exp[1q n (x 1 -x 2 )] (form 2) 

n-1 

2 [|C| 4 /x 2 (N-l)]iN[l + ZJt i |b(k,l) | 2 ] 2 
N/2 ^ * 

- (2/N) £ Cl + 2£j? fJt |b n (k,l)| 2 cosq n (x!-x 2 ) 
n*l 

+.. £ 4 , |b n (kJ)| 2 |b n (lTM , )| 2 cos(q n -q' n )(x 1 -x 2 )]} 

k,l ;k' , 1 ' 

2 [|C| 4 /x 2 (N-1)](N[1 + rf.,|b<r.1)| 2 ]2 -1 -2£iT, 4 |b(k,l)| 2 cosq( Xl -x 2 ) j. (8.9) 

Diis form of the pair-correlation function includes the 1 /N factor which 
multiplies the variable (noise) part. The crucial point in the derivation 
of the 1/N factor was a elimination of the rapidly oscillating terms exp(K«- 
K«)(xi-x*) with K a * K a present in the second form of Eq. (8.9) above, an 
elimination indicated through the approximation sign connecting the second 
form to the third form above. Indeed, since K a differs from K» by much more 
then the momentun change corresponding to the emission of a infraquantun, 
these terms will have a very fast oscillation, and will not yield any low 
frequency noise. Since they are also small in magnitude, they are negli¬ 
gible. This provides the important reduction of the noise term by a factor 
N. In the last form of Eq. (8.8) the first two terms are constant and 
large, and do not yield any rapid oscillations which would justify elimina¬ 
tion of any cross terms. 

In Eq. (8.9) we have used again the form of the single-particle wave 
functions given in Eq. (5.1) or (5.2), with independent sets of phases 







present in the bremsstrahlung energy loss parts of each particle. Had we 
used a random time-shift in each of the single-particle wave functions, i.e. 
just a random initial time constant, or, equivalently, a random space-shift- 
constant, all results would have been exactly the same. In fact, the random 
shift should better describe the initial Poisson distribution of the incom¬ 
ing particles which are scattered. As we have seen in Sec. IV, the random 
phase set in the bremsstrahlung energy loss parts is the set of random ini¬ 
tial phases of the electromagnetic field oscillators. Therefore it should 
come in the same way for all particles. However, the random shift will 
eliminate expectations of the cross products b a *(k,l)b»(k,l) with m ^ n 
just as the sets of random phases used by us did; if P is for instance a 
random space shift, these cross terms yield contributions to the pair- 
correlation function of the form <| b(k,l) | * cosq(xi-x» + p)> = 0, where 
the average is with respect to p. We present this observation here as an af¬ 
terthought, because the sets of random phases generated by a shift in time 
or space will appear to be random, but will still contain some correlations. 
Hie point we make here is that these correlations have no effect on our cal¬ 
culations, so we can continue to use random phases for our purpose. There 
may be, of course, some differences in higher-order correlations which we do 
not consider here. 

In the last form of Eq. (8.9) we have neglected the higher-order 
term. Using Eq. (4.16), we can write the pair-correlation function for fer¬ 
mions in the form 

A(x!.x 2 ) - (|C| 4 /x 2 )tl + [2/(N-l)]Ee,l|t>(ic,l)| 2 [N f t, cosq(xi-x 2 )]}. (8.10) 
dividing again the variable part by the constant term, and neglecting small 
constant terms, we obtain for the fractional spectral density of the fermion 
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current and cross section fluctuations 


Sj(k)dk/j 2 - 2aAdk/k(N-l) - Sj(f)df/j 2 - 2oAdf/f(N-l). 


(8.U) 


IX. DISCUSSION 

As we have seen in Eqs. (2.5a) and (2.5b) the second-quantized cal¬ 
culation of the pair-correlation function for bosons has yielded just the 
autocorrelation of the single-particle probability density which Me had used 
in the initial publications 1 • *. For fermions the same result was obtained 
with a minus sign which corresponds to a 180* phase shift. The latter was 
interpreted as the result of the Pauli exclusion principle which prohibits 
the very close location of two identical fermions of the same spin. 

However, this affects the pair-correlation function only at very small rela¬ 
tive distances of the particles. For the larger distances and lower 
frequencies or wave numbers important for 1/f noise, this difference, which 
appears in the case of fermions, is negligible. We conclude that the simple 
autocorrelation function calculated for the single-particle probability den¬ 
sity in earlier papers on quantum 1/f noise is required by any second- 
quantized calculation which takes into account the quantim exchange between 
identical particles automatically. However, there is nothing special about 
the second quantization beyond this automatism and its neatness; the same 
results are obtained also without it, by properly symmetrizing the wave 

functions, as we have seen directly in Sec. VIII.. Furthermore, quantum l/f 
noise *a* also derived without Che A PS PS WF. 1 * 







At first sight the autocorrelation of the single-particle probability 
density function can not be considered as the expectation value of an ob¬ 
servable quantity in agreement with the principles of quantum mechanics, be¬ 
cause it is of fourth order with respect to the wave function. However, as 
we have seen, this is required by a correct calculation in second quatiza- 
tion. To explain the fourth-order dependence of the pair-correlation func¬ 
tion on the wave function without using the method of second quantization, 
we mention that the nonuniformity of the distribution of the outgoing par¬ 
ticles is described by the two-particle wave function, which in turn is a 
properly symmetrized or antisymmetrized product of two single-particle wave 
functions. Taking the module squared of this two-particle wave function,we 
obtain the probability denr«.'tjr of the relative positions of the par¬ 
ticles, or the pair-correlation function* which therefore will be of fourth- 
order in terms of the single-particle wave functions. This fourth-order de¬ 
pendence which was present in all publications on quantum 1/f noise 1- ®, is 
therefore in agreement with the principles of quantum mechanics, and is jus¬ 
tified by the detailed many-particle calculation, in spite of the 

IS 

acrimonious criticism it has recieved. We have also provided for the first 
time in this paper a formal derivation both of the quantum 1/f cross cor¬ 
relations and of the 1/N factor of the empirical Hooge relations, in perfect 
agreement with the previously introduced 1/N factor 1-5 . 

In conclusion, we have shown how the method of second quantization 
can be applied to derive the pair-correlation function of the conventional 
quantum 1/f effect in space and in time. Both the pair-correlation function 
and the power spectral density are similar to the previously published ex¬ 
pressions based on more elementary derivations. Rather than contradicting 
the existence of quantum 1/f noise, the second quantization method confirms 
the quantum 1/f effect in a brilliant way, and convincingly justifies my 
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previous use of four single-particle wave functions in order to describe the 
self-interference of the particle fields, on the basis of quantum exchange 
effects between identical particles. 

Finally, we mention that the inclusion of the soft photons into the 
final state is not required, if the electromagnetic field modes are not part 
of the system studied, but are considered external, being included as a 
time-dependent external force field in the Hamiltonian. This external field 
description’•■ involves a time-dependent Hamiltonian and a formal lack of 
energy conservation, but is well suited for the derivation of quantum 1/f 
noise and of most other electrophysical and electronics problems, as I have 
pointed out earlier. The present paper shows that the quantum 1/f effect 
can also be derived in the language of second quantization, with the ex¬ 
plicit inclusion of the emitted bremsstrahlung photons into the final state. 

The equivalence of the wave number spectrum with the frequency 
spectrum was proven here by neglecting coherent state quantum 1/f noise, 
which was discussed elsewhere 11 > 11 . Indeed, if the coherent state quantum 
1/f noise is included, fluctuations caused by virtual photons in the final 
state must be included. We recall that coherent quantum 1/f noise is caused 
by the energy uncertainty introduced by the coherent state of the 
electromagnetic field of a physical charged particle. Both coherent and in¬ 
coherent quantum 1/f noise can be derived in the interaction picture. In 
the interaction picture coherent quantum 1/f noise is obtained from diagrams 
which do not involve scattering of the particles, while the conventional (or 
incoherent) quantum 1/f noise considered in this paper arises from diagrams 
of scattering with the simultaneous emission of photons. We have not in¬ 
cluded coherent quantum 1/f noise in our calculation here, because we have 
considered only the noise in the final state long after the scattering 
process, with the coupling to the electromagnetic field disconnected, and 
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with the Heisenberg field operators <J(r,t) in free motion. Without this 
restriction, the quantum 1/f noise spectrum in frequencies would have been 
larger than the wave number spectrum. However, in very small electronic 
devices only conventional quantum 1/f noise seems to be presentM. 
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appendix ft 

THE SINGLE-PARTICLE WAVE FUNCTIONS 

1. Introduction 

The single-particle wavefunctions used in this paper make sense physically 
and were first Introduced in 1975 on physical grounds 1 . The objective of this 
appendix Is to show how they can be derived formally. The single-particle wave 
functions were derived In Sec. IV by treating the electromagnetic field 
classically, and by writing the final result In terms of photon creation and 
annihilation operators. This last step implies the correspondence principle. 

A more rigorous derivation should follow perturbation theory, e.g. the 
calculation of Yennie, Frautschi and Suura 15 reformulated by Chung 16 in order to 
provide a scattering matrix element free of infrared divergences. The 
reformulation by Chung shows that the price we have to pay for a finite matrix 
element Is a more complicated set of incoming and outgoing states describing the 
scattered particles, with the inclusion of a coherent state of the 
electromagnetic field Into the notion of charged particle. This new picture was 
first suggested by Dollard 17 In the language of nonrelatlvlstlc quantum 
mechanics, and was successfully generalized to the relativistic case by Kulish 
and Faddeev 1 * 5 , and by Zwanziger 1 ^. This new picture is intermediary, between 
the usual Heisenberg picture and the usual Interaction picture, because the 
offending (long-range) part of the electromagnetic interaction has now been 
included into the unperturbed hamiltonian. 

Following Chung 16 , we obtain a finite matrix element by writing the initial 
state of a charged particle involved in a scattering process, with initial 
4-momentum p-j and final 4-momentum pf, in the form 

| > 1 - exp(-(l/2)i ktJl |S- t) (k) | 2 } 

ex Pl E k,i, sj^(k)e(*)(k)a k U) + } MPi)> 
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« 11 - (l/2)I k . ll |Sl U) (k)l 2 Jtl ♦ Ik,! s{* ) (k)e(t)(k)« k (»)+j| t (p,)>, (Al.l) 
where the summations over k can be transformed In Intergrals by multiplication 
with (L/2ir) 3 as In Eq. (5.4), and where 

Si (A) (k) - ie/2(2ir) 3 k 0 iP 1 eC*)/kp 1 (A1.2) 

is a function which depends on the momentum pj of the Incoming particle of 
charge e. The state |(p^)> Includes the bare particle and the vacuum on which 
the creation operates aj^ act. Here e(*)(k) is the polarization 4-vector of 
the mode k,i. The physical charged particle defined by the first form of Eq. 
(Al.l) contains the bare particle and its field which Is In a coherent state, 
with an Indefinite number of photons In each mode and therefore with energy and 
momentum which are not sharp, l.e., not well defined. The last form of (Eq. 
Al.l) is only an approximation obtained by expanding the exponentials. The *zrm 
of $j(*)(k) given by Eq. (A1.2) is correct only in the limit k+0 of Interest for 


The final state is obtained by replacing p-j with In Eqs. (Al.l) and 
(A1.2). The photon Is assumed to have a finite mass x which Is allowed to 
approach zero at the end of the calculation. 

It Is easy to verify that if the initial and final particle states are 
chosen as shown In Eqs. (Al.l) and (A1.2), all infrared divergences cancel 
already at the level of the matrix element Mf^, and not only later at the level 
of the process rate or cross section as would be the case had we started with 
the usual notion of particle, without the long-range part of the field included. 
It is easier to do this proof in two steps: first to second order only, in 
order to get acquainted with the formalism; then up to any order. Our treatment 
will be simplified and schematized as much as possible, in order to focus on the 
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Infrared divergences are of two kinds, arising from virtual photons and real 
photons. The former are photons which boomerang back to one of the charged 
particles Involved In the contemplated process. Real photons are emitted (or 
absorbed) for good. Let Mq be the basic matrix element calculated for an 
arbitrary process (e.g., scattering of charged particles on a fixed potential) 
without Including the coupling to photons. In fact, however, the process may 
also happen with any number of virtual photons. Adding up the Infinite series 
of processes (or Feynman diagrams) corresponding to virtual processes of any 
order, one obtains^*® the matrix element 


where we have neglected a small nondlvergent term, and 

, , , f (2(Pf-k) (2(Pi-k)-|(2) 

R.M) * (e2/4V)l( k 2^2)l/2 ^ (2pf-k _ x2) - (2p) . k _ x) J (A1.4) 

- -(l/*)tk,*l4 t) <k)-Si U) (‘)|2 • 

Here the rectangular bracket Is a 4-vector and the exponent In brackets 
designates the scalar product of the four-vector with Itself. For X 3 0 the 
Integral derived from the sum in Eq. (A1.4) shows a logarithmic divergence which 
has been exponentiated In Eq. (A1.3) by summing the the virtual photon processes 
(diagrams) of all orders. The value — of Eq. (A1.4) tells us that M' » 0 in 
the limit X * 0, and therefore there are no processes without the simultaneous 
emission or absorption of any (observable or unobservable) real photons. 

To get a finite result, we multiply the matrix element M 0 which includes all 
virtual photon contributions with the sum of all real photon contributions. To 
second order, this yields 
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M" Ll-(l/2) E |s{ 4) (lc)| 2 jll-(l/2) E |Sf A) (k)| 2 j 

k,i k,i 

U ♦ £ Sp’tkJS^’tkJUl + E sJ l) (k)S (t) (k)Jll - £ s{* , (k)S U lk)J ♦ nM 0 
k,i k,i k,t 

- H 0 [l + n ♦ aB - (1/2) E(|Si | 2 + |Sf | 2 -2sfsf) + E|S k | 2 j 

- MqLI + n + oB + a§J. (A1.5) 

Here In the first form on the r.h.s. we have first Included the normalization 
factors present In the In and out coherent states (Al.l). The following factor 
(third rectangular bracket) Is the scalar product of the In and out coherent 
fields, describing a transition without Interaction, The next factor describes 
the emission of the final coherent state field Sf^(k) In the scattering 
process, and the last factor corresponds to the absorption of the incoming state 
coherent field. Note that the last three rectangular brackets have just been 
taken over and copied from the usual treatment of infrared divergences, just as 
we took over the virtual photon contribution. In the second form on the r.h.s. 
we consequently limited ourselves to the second order In e/(hc)*/ 2 and carried 
out the multiplications. We omitted the upper polarlzatln label and wrote the 
argument k as a superscript, also omitting the indication of the summation 
variables k and i under the sign of the sum. Finally we introduced the 
notations 


Sf°(k) - sj* ; (k) = S(*)(k) aS k ; off - (1/2)e|S k | 2 , (A1.6) 

In order to save space. The term M 0 n stands for a non-infrared-divergent 
contribution. We notice from (A1.4)-(A1.6) that for low k values the summands 
(or Integrands) In aB and aB are Identical, but with opposite signs in the limit 
A-*0. Therefore, the infrared divergences cancel already at the level of the 
matrix element. 













The coherent states of the field of a physical charged particle defined by 
Eq. (Al.l) are vectors in a separable Hilbert space. The result In Eq. (A1.5) 


is not affected by a translation In this space, which amounts to adding or 
subtracting the same vector y k from all states. Indeed, subtracting y£ the 
amplitude S k * Sf - S* remains the same, and we get 

M“fi a M o eoB Li-(l/2)E|Si-y k | 2 jil-(l/2)2:|Sf-y k | 2 j 

Ll-ReE(Si-y k )*(Sf-y k )j[l+ReE(Sry k )S k j[l-ReE(slf-y k )S k j + n M 0 

« M 0 (l+n+aB+aB) (A1.7) 

as before. 

To demonstrate the cancellation of divergences to all orders, we regroup the 
perturbation theory series in the form 

Mfi * M o <0|expLS^a k jexpLES k (a k -a k )JexpL-ES l [a k J |0> 

expL-(l/2)E(|sf| 2 + |Sf I 2 )J + Moh' (A1.8) 

where M 0 n' Is a nondlvergent part. 

Applying the relation e^ + ^ * e P e Q e _ (l/2)CP.Q] ynth p 3 js k a k and Q * P + , we 
obtain 

M f1 » M o e aB <0|expESfa k expEsl[a lc |0><0|expES^a k expES k a) ( |0> 

<0 |exp-cS l! a l( expi:S l JaJ|0> * «„«' ' (AK9) 

The virtual phonon factor e 00 arises from automatically, and has an 

upper Integration limit roughly equal to the energy E of the particle considered 
and involved in the scattering process. The integrals in e a B, however, should 
be divided In a part below the photon observation threshold e 0 ■ hf 0 and a part 
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above, going also up to E, which contains the part of the matrix element 
occurlng with bremsstrahlung. The parts below c 0 cancel each other In «B and In 
oB, demonstrating the disappearance of all divergences. Bringing the rest of 
e® 8 to the denominator and neglecting the difference between B and §, we can 



The sum with a prime is over |k| > e 0 only, l.e., is restricted to observable 
photons. Here we have separated in the numerator the unity which represents the 
non-bremsstrahlung part of the rate, as we did earlier* to second order. 



3. Single-particle wave functions 

The discussion and derivation outlined above suggests Eq. (Al.l) as the 
correct choice for the incoming and outgoing physical charged particle states. 

For the out state the subscript i is replaced by f. However, conventional 
quantum 1/f noise (effect) Is defined as the cross section (or process rate) 
fluctuation arising due to the small bremsstrahlung effects associated with the 
process considered. The current fluctuations caused by the general uncertainty 
In the energy and momentum of the physical charged particle states have always 
been called coherent 11 * 1 ^ state quantum 1/f noise (effect), and will be studied 
separately in another paper****^. Since we are interested only in the noise 

contribution introduced by the collision at hand, we perform a translation in 

k 

the space of principal (coherent) state vectors, by taking y * in Eq. (A1.7). 

The Invariance of matrix elements with respect to any translation was verified 
above In second order, and is manifest in Eqs. (A1.10) and (A1.9) to any order, 

f k 

because S k is an invariant difference of two vectors. Subtracting S-j from the lkl>£c 
amplitudes of all vectors we obtain a bare particle incoming state with no obsarvaile 
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photons (|k|>e ), and an out state 


|> f = exp{-(l/2)E e>|l |S l (k)|2) 


exp(£ Eft s(*)(k)e( t )(k)e- i 5' I - f )a|L‘) t ) | ?( p ( )> 
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tl-(l/2)r k - t JS^k)][l + E j j jl S ( ^(k)e' i P ( ^ ) a^ )+ ]|4-(p f )> (A1<11) f 


Here we have introduced also the recoil phase factors e "'' c l( x ‘ r ) which describe 
the momentum loss q * kc/v of the outgoing paricles of momentum p f = mv when 
a photon of momentum k is emitted. These recoil phase factors are neglected in 
the independent boson model which ignores the changes caused by the emission 
of individual photons; we have neglected these factors so far, but are restoring 
them here. The particle-specific constant r has to be introduced because the 
single-particle states |> f are not exact momentum eigenstates, and are there¬ 
fore not exactly invariant with respect to translations. Individual single 
particle wave functions may differ by such a random translation in space or 
time regardless of whether they are localized or completely unlocalized in 
space and time. For momentum and energy eigenfunctions, such translations 

introduce only a general phase factor which drops out of any physically relevant 

In Cj/. (Al.U) JTJatt n c r drep o ut. J J thifi T i/jradial COtrJol ate X i n the te^u/wak/ltfe. , 

calculationVFor a spherical wave the translation is replaced byY"a time-shift, 

or time delay; if this requires more than a general phase factor, a randomrad.'a/(cr 

timelshift must be introduced for each scattered particle in the exponents. 

With the notation S^(k)e^e lc,r = b(it,z), Eq.(Al.ll) becomes similar to 

Eq.(5.2) which was used in the main text of this paper. The main difference 

is the normalization factor in (Al.ll), which was suppressed in Eq.(5.2). The 

phase factor e iqr (or e^ r for spherical symmetry), with constant, particle- 

specific r, present in b(l<,*), is the random phase factor introduced by 
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us in the main text for each particle, which led to the 1/N factor in 
Sec. VIII. This derivation of the single-particle wave functions 
treats the electromagnetic field quantum-mechanically. while the deri¬ 
vation in Sec. IV was semiclassical. 

In the present paper we have considered only one isolated scatte¬ 
ring event. We also note that the subtraction of S- 1 changes the noise 
given by the final state. The corresponding separation in conventio¬ 
nal (incoherent) and coherent state quantum 1/f noise seems to be con 
venient both for conceptual or didactical, and for practical or calcu- 
lational reasons, connected with the method of compounding the noise 
contributions from successive scattering events. As was indicated 
earlier 1-0 , this addition of contributions from many successive, closely 
spaced, scattering events gives a result close to the contribution, 
calculated here, from a single representative scattering event with a 
correction factor which is usually close to unity. 








APPLICATION OF THE QUANTUM 1 /i THEORY TO MIS INFRARED DETECTORS 







I. INTRODUCTION 


The utilisation of infrared detectors at larger wavelengths and in the staring 
mode has made the limitations imposed by i/f noise on the performance of MIS 
structures more stringent and conspicuous. On the other hand, the development 
of the quantum i/f theory has provided us for the first time with a possibility to 
predict, model and calculate from first principles the various i/f noise 
contributions affecting the many currents and processes which are of importance 
in the operation of infrared detectors. 

The mam purpose of the present report is to apply this new knowledge of i/f 
noise to MIS structures working as infrared detectors. We start here with a brief 
general description of quantum i/f noise and continue in Sec. II with the 
inventory of various components present in the current of MIS infrared detectors. 
In Sec. Ill we analyse the quantum i/f noise associated with the currents and 
processes discussed in Sec. II. In Sec. IV we compare the i/f noise components in 
magnitude and determine their impact on the performance of MIS detectors. 
Finally, in Sec. V we discuss the resulting i/f noise limitations, and point out 
some possibilities of reducing the quantum 1/f noise. Those familiar both with 
quantum i/f noise and with infrared detectors may now go directly to Sec. III. 

Quantum 1/f noise*”^ * s a fundamental fluctuation of physical cross sections 
and process rates, caused by the infrared - divergent coupling of current carriers 
to low frequency photons and other infraquanta. The physical origin of quantum 
1/f noise is easy to understand. Consider for example Coulomb scattering of 
electrons on a center of force. The scattered electrons reaching a detector at a 
given angle away from the direction of the incident beam are described by 
DeBroglie waves of a frequency corresponding to their energy. However, some of 
the elecrons have lost energy in the scattering process, due to the emission of 
Bremsstrahlung. Therefore, part of the outgoing DeBroglie waves is shifted to 









slightly lower -frequencies and inter-feres with the main, non-Bremsstrahlung, 
component, yielding beats. These beats present in the probability density along 
the direction o-f the scattered beam will be noticed in the detector as low 
frequency current fluctuations, and will be interpreted as fundamental cross 
section fluctuations. Although the wave function v |' of each carrier is split into a 
Bremsstrahlung part and a non-Bremsstrahlung part, no quantum 1/f noise can be 
observed from a single carrier. A single carrier will only provide a pulse in the 
detector. Many carriers are needed to produce the i/f noise effect, just as in the 
case of electron diffraction patterns. While incoming carriers may have been 
Poisson distributed, the scattered beam will exhibit super - Poissonian statistics 
or bunching due to quantum i/f noise. The quantum 1/f effect is thus a two - 
particle effect, best described through the two - particle wave function and two - 
particle correlation function. 

Let us estimate the magnitude of the quantum i/f effect by starting with the 
classical (Larmor) formula qv 2 /3c 3 for the power radiated by a particle of charge 
q and acceleration"^. The acceleration can be approximated by a delta function 
^t) =^"ve/'<t) whose Fourier transform^"v is constant. The one - sided spectral 
density of the emitted Bremsstrahlung power 2qCa,"v) 2 /3c 3 is therefore also 
constant. The number 2q(p!v> 2 /3hfc 3 of emitted photons per unit frequency 
interval is obtained by dividing with the energy hf of one photon. The probability 
amplitude of photon emission C2q&v) 2 /3hfc 3 ] i/2 is given by the square root of 
this photon number spectrum, including also a phase.factor. The beat term in the 
probability density l^l 2 is linear both in this Bremsstrahlung amplitude and in 
the non - Bremsstrahlung amplitude. Its spectral density will therefore be given 
by the product of the squared probability amplitude of photon emission with the 
squared non - Bremsstrahlung amplitude which is independent of f. The 
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resulting specti al uensity of -fractional probability density fluctuations is 
obtained by dividing with I'f'l 4 and is therefore 

lfl“ 4 SI l 2 (f) = Sq 2 (Xv) 2 /3hfNc 3 = QA/iU = r 2 Sj(f), (1.1) 

where ot = 2e 2 /hc = 1/137 is the fine structure constant ando^A = 4q 2 k3'^)“/3hc 3 is 
Known as the infrared exponent in quantum field theory, and is known as the 
quantum 1/f noise coefficient, or Hooge constant, in electrophysics. 

The spectral density of current density fluctuations is obtained by multiplying 
the probability density fluctuation spectrum with the squared velocity of the 
outgoing particles. When we calculate the spectral density of fractional 
fluctuations, the velocity simplifies and therefore Eq. (1.1) also gives the 
fractional spectrum of current fluctuations Sj(f>, as indicated above. The quantum 
1/f noise contribution of each carrier is independent, and therefore the quantum 
1/f noise from N carriers is N times larger; however, the current j will also be N 
times larger, and therefore in Eq. (1.1) a factor N was included in the denominator 
for the case in which the cross section fluctuation is observed on N carriers 
simultaneously. 

The fundamental fluctuations of cross sections and process rates are 
reflected in various Kinetic coefficients, such as the mDbility/^fand the diffusion 
constant D, the surface and bulK recombination speeds s and recombination times 
Zt the rate of tunneling and the thermal diffusivity. Therefore, the spectral 
density of fractional fluctuations in all these coefficients is given also by Eq. 

(1). This is true in spite of the fact that each carrier will undergo many 
consecutive scattering processes in the diffusion process. The quantum i/f noise 
in the mobility and in the diffusion coefficient is practically the same as the 
quantum 1/f noise in a single representative scattering event which limits the 
mobility or the diffusion coefficient. 
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Due to the rapid relaxation of concentration fluctuations, the quantum 1/f 
fluctuations of scattering cross sections will only be reflected b'/ the 
fluctuations of the mobility and the diffusion constant of the carriers, and not by- 
fluctuations in the concentration of carriers. 

For large devices the concept of coherent state quantum 1/f noise was 
introduced 1 12 ; see j ,-, kzftor details. In this case the Hooge parameter 
may be written 

* H * <*H>coh a 2 a 4T S 4.4 10” 3 (1.2) 

where ctf = 1/(137) is the fine structure constant. This is of the same order of 
magnitude as the empirical value <Xjj = 2 10“ 3 -that Hooge found for long devices. 
It is therefore proposed that Hooge's emprical value for°^j^ is due to coherent 
state quantum 1/f noise, so that it has a very fundamental origin. 

For small devices (e.g.» of size L < 10/<m) we apply conventional, or 
incoherent 2 ”®, quantum 1/f noise which is just the cross section fluctuation 
introduced above in Eq. (1.1). In that casec^ may be written 

p<H s ^H’incoh = <4<*/37r>CUv) 2 /<c 2 )D, (1.3) 

where av is the change in the velocity of the carriers in the interaction process 
considered. This expression holds for any 1/f noise source describable by 
fluctuating cross sections. Since usually (av* 1 /c 2 )<3fl, except for carriers with a 
very small effective mass, we now have o(j^ < 3.1 10”®. This may explain the low 
values of (in the range of 0 ^ = 10”^ - 10“^) for very small devices. In 
between one can introduce a parameter s = f(L/L 0 ) where L D is a characteristic 
size and write 12 . 

*H * ( Vincoh C1/<1 + 5,3 + < V hWs /(1 + *»• 


(1.4) 
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w:th s« i -for L/L C <?1 and s^>l -for L/L D >)i. Aczording to this rough 
approximation 12 , L r , = lOy^m -for samples with a concentration c of carriers of 
10^cm - ^ an( j vanes proportional to c _1//t: . This describes the transition -from 
Eq. (1.2) to Eq. (1.3) when one goes to devices with smaller and smaller sizes, as 
we shall see below in Eq. (3.21). 

When we apply Eq. (1.1) to a certain device we -first need to find out which are 
the cross sections which limit the current, and then we have to determine both 
the velocity change &v of the scattered carriers and the number N of carriers 
simultaneously used to test each of these cross sections. The next section deals 
with the first problem, i.e.» is dedicated to the study of currents in MIS infrared 
detectors. 


II. CURRENTS IN MIS DETECTOR STRUCTURES 

MIS detectors 9 are different from photovoltaic detectors, because they do not 
contain a pn junction obtained by inhomogeneous doping, and use an insulated 
field plate, or gate, placed on top of a homogeneously doped narrow - bandgap 
semiconductor. The gate is used to control the surface potential, driving the 
semiconductor surface into deep inversion. The field of the induced quasi - pn 
junction obtained under the surface of the gate in the homogeneous semiconductor 
material is used to separate the carriers generated by photo-electrically induced 
band to band transitions just as in a photovoltaic device. The MIS infrared 
detector is therefore similar to a capacitively coupled photovoltaic detector, 
without the inconvenience of inhomogeneous doping processes. 

MIS detectors are operated in the pulsed regime by applying the gate potential 
which creates the inversion under the surface for a finite time only, and by 



applying subsequently a potential which flattens the energy bands r 


surface and releases the carriers which had accumulated from photo electric 




I 
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effect and dark current processes during the preceeding interval. The electrical 
signal obtained when the carriers are released, i.e., during readout, is 
prcooruonal to the number ot' carriers accumulated, and therefore to the total 
current supplying the inverted volume under the surface with minority carriers 
from the bulk and from various thermal and photo electric processes in the 
depletion, inversion and surface regions. This electrical signal is used in order 
to determine the flu>: of infrared radiation. For this determination, however, the 
dark current contribution needs to be subtracted first. 

The dark current is the current supplying the potential well, defined by the 
inversion region under the surface, with minority carriers in the absence of the 
applied infrared flux. Any low - frequency fluctuation in the dark current will be 
interpreted as a fluctuation in the major infrared flux signal. Therefore 
fluctuations of the dark current at frequencies below the readout frequency will 
limit the performance of infrared detectors. In the pulsed mode of operation 
considered here, the dark current is monitored only during the inverted phase, 
when carriers are accumulated in the potential well. Therefore the cross sections 
and process rates which control the intensity of the dark current are not observed 
continuously either. Nevertheless, the quantum 1/f fluctuations of these cross 
sections and process rates will be the same as if we would have observed them 
continuously. Indeed, the changes in the incoming flux of electrons testing all 
cross sections and process rates in the semiconductor is only slightly affected by 
the applied gate voltages, and is present also in thermal equilibrium. This 
independence of 1/f noise on the continuous or discontinuous character of any 
applied bias has been experimentally verified 16 during the last 2 decades, and 
has been found to be in agreement with the interpretation of 1/f noise in terms 
of fundamental resistance fluctuations. Although the experimental verification 
was performed cn fluctuations in conduction only, from the concept of Quantum 
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1 /-f noise we know that the similarity of quantum l/'f noise in the continuous and 
pulsed regimes should be also tnje tor quantum i/+ -fluctuations in recombination 
cross sections and tunneling rates. 

The dark current has to be subtracted from the total current in a (HgCd)Te MIS 
device to yield the photocurrent. Therefore, the minority - carrier dark current is 
the single most important parameter tor the operation o-f MIS devices as 
detectors infrared radiation 9 . This applies both to operation of MIS devices m 
the thermal equilibrium mode, in which the dark current determines the MIS diode 
impedance, and to operation in the dynamic, or integrating mode, in which the gate 
voltage is pulsed, and in which the minority - carrier dark current determines the 
storage time of the device. The main component of the dark current in narrow 
bandgap HgCdTe is the tunneling current via bandgap states 9 , which can also be 
considered as an electric breakdown effect. In general, the tunneling current 
occurs both through band to band transitions and through intermediary states. 

The band to band tunneling current through a simple triangular barrier is 

J tb = (q 3 E < ^/4X 3 h 2 )(2m*/E g ) 1/2 expC4(2m*)E g 3/2 /3qhE3, (2.1) 

where E is the electric field associated with the barrier, and E g is the bandgap. 
The electric field can be approximated by the electric field at the semiconductor 
surface 

E 4 .= (2qn 0 (^/££ 0 ) 1/2 , (2.2) 

where represents the empty well surface potential, and n Q is the doping 
concentration. Substituting this value into Eq. (2), with m*/m c , = 7 10 -2 E C , we 
obtain 9 

J tb = 10" 2 n o 1/2 (f) s 3/2 e>:pC- 4.3 10 10 E g 2 /(n o <fe) 1/2 : A/cm 2 , (2.3) 





m\ 













where n 0 is in cm -3 and Eg,1^ in volts. Therefore the tunneling current is 
strongly dependent on the bandgap and also depends on the doping concentration 
and the surface potential. 

Experimental values of the tunneling current are usually larger than Eq. (2.3) 
because of the additional effect of tunneling via bandgap states. This effect is 
particularly important in n - type des'ices 1 ®. Indeed, in n - type devices the 
applied gate voltage is negative in order to produce depletion at the surface. The 
energy bands are therefore curved upwards at the surface, and transitions of 
electrons from the valence band to Shockley - Read (SR) states at the middle of 
the bandgap, as well as the subsequent transitions from these states to the 
conduction band are facilitated by the presence of many defects right at the 
surface of the semiconductor. In p - type devices the similar indirect tunneling 
processes occur father away from the surface, because in this case the bands are 
curved downwards at the surface, and transitions of electrons from the valence 
band to the centers at the middle of the bandgap, as well as the transitions from 
the centers to the conduction band well at the surface, occur right where the 
curvature begins, i.e., further away from the surface. We conclude that in p - 
type devices there will be fewer SR centers active in indirect tunneling, and 
therefore the tunneling current J c via SR centers at a given teperature and a 
given applied gate voltage will be smaller. The tunneling current will be further 
reduced in p - type devices due to the lower density of states present in the 
surface potential well due to quantisation of the motion of the electrons in the 
potential well at the surface. The reduced values of the dark current in p - type 
devices correspond to higher values of the breakdown field in these devices. The 
best measured valued of the breakdown field in 10cutoff p - type devices is 
in e>;ess of 1.0 V/^m, whereas that for n - type material of similar bulk defect 
quality is 0.5 V/um. On the other hand, the minority carriers diffusion current is 
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larger in p - type devices due to the smaller mass and higher diffusion constant 
of electrons compared to holes. The advantage of p - type devices is therefore 
considerable only in the case of very narrow bandaap and very long cutoff 
wavelengths. We shall therefore consider both the case o-f p - type and n - type 
devices. The large diffusion current present in p - type devices corresponds to 
the large value of the diffusion length of electrons and can be reduced by thinning 
the devicei i.e., by reducing its thickness well below the diffusion length. 

In general the dark current can be written in the form 

J d = J diff + J deo + J s + J tb + J tc + J tsc + J b + {2A) 

The eight terms on the right hand side correspond to minority carrier diffusion 
from the bulk, generation from SR centers in the depletion region, generation 
from SR centers at the surface, band to band tunneling, tunneling via SR centers, 
tunneling via surface centers, recombination on the back surface and photoelectric 
generation by the thermal radiation background flux^. We shall give the 
formulae 10 which determine these terms below, also including an e>:ample of their 
calculation in a p - type device 

J dif = <qn i 2 /p 0 )Ck’I/^qT a ] 1/2 = (1.6 10~ 19 C 36 10 24 cm“ 6 /10 15 cm~ 3 )C0.01 V 
1.5 10 5 (cm 2 /Vs)/10~ 6 s: 1/2 = 2 10“ 4 A/cm 2 , (2.5) 

J depl = P^W/2To = C(1.6 10“ 19 C 6 10 l2 cm -3 2 10 _4 cm /(2 10“ 6 s)D 
= 10 _4 A/cm 2 (2.6) 

j c - J Q = qnjS/2 = [<1.6 10 _19 C 6 10 12 cm -3 20cm/sD = 1.3 10 _5 A/cm 2 , (2.7) 

J tc = 10 _13 (N r 4/E g ) exp[-5.3 10 6 E g 2 /E3 = 10 -i3 <4 10 12 cm _3 0.2V/0.062V) 
e>:pC-5.3 10 6 <0.062V> : /3000 V/cmD = 8 10“ 4 A/cm 2 , (2.8) 
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and Jtb was S lven by Eq* (2.3), yielding 10“ 9 A/cm 2 . Here n, = 6 lO^'cm" 3 is the 
intrinsic carrier concentration, and the concentration of holes was caken to be p D 
s 10 15 :m“ 3 . The mobiiity/Oi = i.5 10 i5 cm 2 /Vs as well as the life time X A - i0“ 6 s 
or the minority carriers have to be replaced by,^p and2y> tor the case of n - type 
devices. The surt'ace recombination speed was denoted by s = 40cm/s, and the 
concentration or intermediary states effective -for tunneling was denoted by N p * 
A iO* 2 cm“ 3 . The bandgap considered was E_ = 0.065eV, the surface Dotentiaifs * 
0.2V and the electric -field below the surface E = 3000V/cm. All numerical values 
have Deen included only as an example and are not characteristic of a particular 
device. The numerical -factors in Eq. (2.S) correspond to p - type HgCdTe with a 
20^ m cutoff wavelength and were taken from the paper of Kinch and Beck 1 ^. For 
n - type devices we also need to include tunneling via surface states. If the 
density of fast surface states is denoted by N +s = (0 12 cm^V 1 , the current 
generated from tunneling via a uniform density of fast surface states across the 
bandgap will be given by 9 

J tsc = -qy5iN fs kT/2q = i.75 10‘ 4 A/cm 2 . (2.9) 

The sum of the first seven currents on the right hand side of Eq. (2.4) must be 
smaller than the eighth which corresponds to the thermal radiation background 
flux, for background limited (BLIF) operation. Although not all terms in the dark 
current are of importance, we still retain them at this point, because their 
quantum 1/f noise may be quite significant, even if the corresponding current is 
negligible. We shall now proceed with the calculation of quantum 1/f noise 
contributions from all these currents. 
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III. QUANTUM 1/f NOISE SOURCES 
1. 1/f Noise in the Diffusion Current 

The diffusion limited dark current exhibit 1/f noise due to 

conventional quantum i/f fluctuations in the scattering cross sections of the 
carriers due to phonons and impurities. We apply the fundamental formula civen 
by Eq. <i.i> for an individual scattering process in which the velocity change is 
given by the thermal energy of the carriers, with the assumption that th? 
collisions are perfectly randomising collisions. If the velocity"^ is rotated by an 
angle $ in an elastic collision, the velocity change isl^js 2v sin(<?/2). Averaging 
over all angles and velocities, we obtain 

/v 2 = 4v 2 sin 2 <6/2) = 2v 2 . (3.1) 

and therefore from Eq. (1.1) we get in thermal equilibrium at the temperature T 
the 1/f noise coefficient 

o^H = <4ot/370<6kT/m*c 2 ), (3.2) 

where we have assumed a Maxwell distribution of velocities. For Hg^Cd^Te 
with :< * 0.2 we have m n * ■ 0.008m o and for x * 0.3 we have m n * = 0.02m o . 
Therefore we obtained s 2 10~^ in the first case and^g = 7.5 10 -8 in the second 
case. 

For the case of Umklapp scattering, which occurs in semiconductors only to a 
limited extent due to the relatively small number of high momentum phonons 
available at the temperature T, the momentum change of the electron is given by 
the smallest reciprocal lattice vector, and therefore^ v = h/am*. We therefore 
obtain the quantum 1/f noise coefficient 


0(. H = <4<*/370<h/m*ac) : 


(3.3) 
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which is much larger than Eq. (3.2), but has to be multiplied with a negative 
exponential which describes the scarcity of phonons with momentum of the order 
of a reciprocal lattice vector. [The negative exponential e”^ /T could be included 
in the current weight ■factor which will be de-fined below in Eq. (43), but we prefer 
to include it here already]. Combining Eqs. (3.2) and (3.3), we obtain -for 
conventional 1 /f noise in the mobility and oi-f-fusion coefficients 

<* H = (4«/3iT)[(6kT/m*c 2 ) + (h/m # ac) 2 exp(-^/T)3, (3.4) 

where & is about halt the Debye temperature -for simple metals, but may be 
higher, or the order o-f the Debye temperature, -for semiconductors. 

The quantum i/f noise considered so -far is known as conventional quantum 1 /f 
noise, and a-f-fects cross sections and process rates. In sufficiently large 
semiconductors samples we expect a larger -form of quantum 1/f noise, described 
in Appendix A and known as coherent state quantum i/f noise. For this type the 
1/f noise coefficient is given by 

cx coh = 2 «W = 4.6 10“ 3 . (3.5) 

The values of the quantum 1/f noise coefficient given by Eqs. (3.1) - (3.5) can 
be used to calculate the quantum 1/f noise which affects the various currents 
listed in Eq. (2.4). We first consider the case of the dark diffusion current of 
electrons from the bulk through the surface barrier in a p -type MIS device, 
similar to diffusion in a n + p junction, because in both cases the current is 
determined by the diffusion of electrons which are minority carriers, against the 
built - in field of a Boltzmann potential barrier into the surface well, and by the 
thermal generation of carriers there. We start with the derivation of the 
mobility fluctuation part of quantum 1/f noise in a n + -p diode. For the MIS 
barrier, just as for a diffusion limitted n + -p junction, the current is controlled by 


n 
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dilution of electrons into the p - region over a distance of the order of the 
diffusion length L = (D n ‘T n ) 1 '' 2 which is usually shorter than the length w p of the 
p - region. If Nix) is the number of electrons per unit length and D n their 
diffusion constant, the electron current at x is 

l nd = - eD n dN/dx, (3.6) 

where we have assumed a planar junction and taken the origin x = 0 in the junction 
plane. Diffusion constant fluctuations, given by K7/e times the mobility 
fluctuations, will lead to local current fluctuations in the intervali x 

Jl nd (x,t) = I_jj{4P r| (>! ,t>/D n . (3.7) 

The normalized weight with which these local fluctuations representative of the 
intervals x contribute to the total current l d through the diode at x = 0 is 
determined by the appropriate Green function and can be shown tD be 
(Ax/L)exp(-x/L) for w p /L » i. Therefore the contribution of the section^ x is 

ckl d <>:,t) = (Px/L)exp(-x/L)I nd <fD n (>:,t)/D n , (3.3) 

with the spectral density 

S^id^.f) = (/lx/L) 2 exp(-2x/L) I nd 2 S Dn (x,f)/D n 2 . (3.9) 

For mobility and diffusion fluctuations the fractional spectral density is given by 
o^Hnd /,f ^ >: ' wher ^-Hnd is determined from quantum i/f theory according to Eos. 
(3.1) -(3.5). With Eq. (3.6) we obtain then 

S^id^x.f) = &x/L 2 ) exp(-2x/L) ieD n dN/dx) 2 o( Hnd /fN. <3,10> 

The electrons are distributed according to the solution of the diffusion equation, 
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=^ hpr /fAN 
since jilu!) = /IN/Ax. 

IX . V) = A I- 


Hnr^ N3 
: 2 /tt. 2 NO:)}Ax 


where ll'ix) = CMtt 


- N_3exp(-x/U and Nix) * !!'(>:) + M_ as be tore. 


I* is easily snown -that -the -fluctuating current</AT(><, t) at the junction is 

©fel'.x, t) =Ja Z r (x, t)exp(-x/L> <3.13) 

so 

Sjit) =0 ( Hnr ie 2 N p L/-f2' n 2 ) J <CN'(x)3 2 /NpM(x)}expi-2x/L)d<x/L) = 
oL Hnr (e V*V f E* 2 u 3 /Uu + l)Ddu =c(, 4nn Cei/-fr r( j[F(a) - F(ae~ w/L )3 ‘3.19) 

where F<a), a, and u have the same meaning as before. 

We can use the similarity of the quantum i/f noise results for diffusion 
current fluctuations caused by mobility fluctuations and by recombination speed 
fluctuations in order to combine both into a single formula 

Sj(f) =(d Hd +^ Hnr )Cel/f£ n 3CF(a) - F<ae~ w/L >3. (3.20) 

In the limit of very short devices (w p « L) the last factor becomes aw/C(a + 1)Lj, 
and in the limit of long MIS devices <w p >> L) it simply becomes F(a). In addition 
we have a current noise contribution s Jb from the quantum i/f fluctuation of the 
recombination speed s on the back surface. 

Sc far we have considered only conventional quantum i/f noise which is 
applicable to sufficiently small devices. In general, however, we must 
interpolate between conventional and coherent quantum i/f noise, according to 
the relation 




‘ C i / < 1 * s) j CZ A/il + Cs/(1 + s)2Z2c£//ii J, 
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where s = E m /E^ = Se^N'/m*^, and N' = nS is the number of the carriers per unit 
length c f the device, s represents the ratio between the magnetic energy per unit 
length and the kinetic energy per unit length ot the device. The quantity e 2 /m*c 2 

* r 0 m/m* can be calculated in terms of the classical radius or the electron r c = 
2.34 10~ 13 cm. Then we obtain s = 2 r.-N'm/m*, i.e., the parameter s represents 
twice the number o-f carriers present in a length o-f the device equal to the 
classical radius ot the electron. We must compare s with71 A, and if s <<7 Ta we 
apply conventional quantu.m i/'f noise, whereas for s >>n A we have to apply 
coherent quantum i/f noise. In general the approximate formula of Eq. (3.21) 
must be used for the transition region. 

The dimensionless parameter s is easy to calculate in any practical case. For 
instance in the case of a MIS device of area 50yU.m x 50/i.m with a concentration of 
carriers of iO^ cm -3 we obtain N' = 2.5 10^/cm, and with m/m* = 50 we obtain s 

* 7. On the other hand, we can estimate A for conventional quantum i/f noise and 
we will certainly find A « 1, because the velocity change of the carriers must be 
much smaller than the speed of light. Therefore, in this case we must apply 
coherent quantum i/f noise, because s >> 7TA. Consequently, in Eq. (3.20) we must 
set 

^Hd + ^Hnr =o6 coh = 4 - 610 " 3 (3 * 22> 

The coherent state quantum i/f noise coefficient thus replaces the total 
conventional Hooge parameter. 

2. i/f Noise of the Recombination Current Generated in the Depletion Region 

The quantum i/f noise of the recombination current thermally generated in the 
depletion region arises from quantum i/f fluctuations of the bulk recombination 













rates :n the depletion region. The difference between the recombination rate R 
and the generation G is given by 


K - C* = Con - n i i ]/C(n + nj)r p0 + (p + Pi)r nQ ] (3.23) 

where nj and pj are electron and hole densities when the Fermi level lies, at ths 
trap level. 1f the trap level lies at the intrinsic level, ,Oj = pj = n^. Moreover, 
c po andT no are time constants tor electrons and holes. If A is the cross - 
sectional area of the junction, the current is 


I ■ Je(R - G)A dx = e f Cpn - n 2 *: 

o J 


:/C(n + n,)r po + <? + n^tno^A d>: (3.24) 


where w is the width of the space - charge region and the trap level is assumed 
to lie at the intrinstic level. 

We now turn to the g - r noise. The time constants Tr po and T nD fluctuate in 
1/f fashion and this produces the quantum i/f contribution to g - r noise. We nc 


where C n and C p are the generation (or combination) rates for a single electron 
and for a single hole, respectively. Consequently 


*Wrno=-<JC n /C n ) 


We now aoply this knowledge to Eq, (3.23) and observe that 


<5(R - G) = C R (:<) - G(>:) 3 

CC(n +n i )r DD (<rC D /C ;? )3 + C<p + Pl >r nD tf'C n /C n >3>/C(n *n t )T + ip +n i )7 nc 3. (3.2?) 
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& I = e j <fCR(x> - G(x)JAdx, 


the noise is 

v,-u‘ 

5 T (-f) = e" J j lR( x) - G«x)3AC(R(>:') - G(x')3A 


CC(n + n i ) 2 r p0 2 5 Cp (x, f)/C p 2 + (p + n^T^S q^a, x'» 

/C(n + n^Tpjj + ip + n i )^ n0 D 2 >dxd:< / , <3.29) 

since 5 C p and 6 C n are independent. 

We now observe that 

pn - n 2 i = n 2 i Ce>:p(eV/2KT) - 1 DCexpteV/2KT) + ID (3.30) 

and that the integrand in Eq. (3.29) has an appreciable value only i-f p~ n - n i 
exp*eV/2KT). By substituting n + ^ = p + r.j = n i Ce>:p(eV/2k F 0 + 13 we de-fine an 

effective width such that 
40- 

eA J t(pn - n 2 i )dx3/C(n + n^Z^ + <p + 

= eA{(pn - n 2 ^)/n^Ce>:p(eV/2KT) + 13}Cw e ^/<r po +i- no )3 (3.31) 

We may thus write 

I s lg r = eAw e ^n i Ce>:p(eV/2kT) - 13/(Z po + 2T no ) 

= teN e ^/(7 po +r no >3tanh eV/2KT (3.32) 

where N pff = AWp^CexpteV/2KT) + 13 is the effective number of hole - electron 
pairs taking part in the conduction and noise processes. This equation is exact 
but not very useful since it contains the unknown parameter w^. 

We now turn to Eq. (3.29) and observe that 
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S Cp («, x', f)/C 2 p = Hp /f)/CR(>:') + GO:')3(Z* po +r no )A}J(x' - >:) (3.33) 


S Cn (>:, x', f)/C 2 n = Cfct Hn /f)/CR(>:') * G(>:'):(T po +T no )A3<f(:<' - >:) (3.34) 

The -fac-tor (r n0 +2T pcj ) enters in because S Cp 0:, fl/C^p and S Cn <>:, y.', f)/C^ n 
must be independent o-f 2"p 0 and 4" no if p~n-nj e>:p(eV/2KT). This yields ( if we 

integrate over the $ function 
V 

Sj(f) = e 2 f <CR(x) - G(:<)j 2 A/(r po +r n o ,CR( " ) + G<>!,3> 

_ tr „ 

<C(n + nj)^p 0 ^ Hp " + <P + "i^no^Hn" 3 '"" + n i ,r po + <P + "i^no^Jdx (3.35) 

We now observe that the second factor in So. (3.35) is practically a constant 
as long as p and n are comparable. We may thus bring that factor outside the 
integral sign and write 


C(n ♦ n^ZTpc&Hp" + <P + "i^no^Hn" 3 
/C(n + ni >r po + <p + n^rnol 2 =<* H /f , <3.36) 

whereC^ ^ is given by 

d~ H = ^po^po + L no^ Hp + <^no /4 po + ^ no^ Hn < 3 * 37) 

We thus have 

'Ut 

SjW = 1 ^ He /f<t no + ~po )3 e J CR<>:) “ G<«J3 2 A/[H(x) + G(x))dx 
*& H e V*<~no + T po )3tanh[eV/2KTJ (3.3S) 

We can now prove Eq.(3.38) in a Hooqe - type formulation Cl D. 
Here we put 

SiW/l^r H^^eff (3.3?) 


M 
sia 
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(3.3?) 










But according to Eq.(3.32>» I or = (eN e ^/Z/tanh eV/2kT, so that 


Sj(-f) =oi H CeI gr /*r)tar,h eV/2kT, T= (r no +2r po ) (3.40) 

in agreement with Eq.(3.38). 

We can also prove Eq.(3.38) -from the following consideration. We write I cr = 
(No^^^l) tanh eV/2kT, where h\ = e/TTfluctuates. We then have 

S^jW/A I 2 =d H /f or S Al (f) =0C H /f(e/ri4l (3.41) 

so thati since the hole - electron pairs are independent 
S : (f> = N eH S 4! (f) = <eN eff dI/fZ# H = C(eI gr /fTbC H Jtanh ev/2kT (3.42) 

whereof is given by Eq.<3.37) and T -XZ no +Z Dc< ). 

The last two approaches are easily extended to other cases; the method works 
as long as a time constant X and an N e ^ can be defined. 

We finally evaluated. ^j p andd H n * r o m quantum 1/f noise considerations C23» 
C3 3. 

d Hn = <4d/3KXAv 2 n /c 2 ) = 4d/37TU2ea(V dif - V) + 3kT3/m* n c 2 } (3.43) 

d H p = (4oi/370(Zv^p/c 2 ) = 4d/37HC2e(l - a)(V dif - V) + 3kTJ/m* p c 2 l (3.44) 
and as a consequence (see Eq.(3.37) and 
^ H = <4d/3TO{C2e<V dif - V) + 6kT3/C(m* n ) 1/2 + (m* p ) 1/2 3 2 c 2 (3.45) 

The problem has hereby been solved. Note that in Hg j_ >: CD >; Te with :: = 0.3! m n 
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3. Noise in the Bur-face Recombination Current 

The surface recombination currnet is a dark current component originating 
■from surface states at the interface between the surface passivation layer and 
the bulk. The recombination cross sections of these states will exhibit quantum 
1/f noise. The quantum 1/f noise coefficient^ ^ reflects the velocity change of 
carriers involved in this process 

oi H = <40t/37t)<2/m*c 2 )C3kT/2 + eU/2 + G.leVl, (3.46) 

where U is the surface potential jump present between the surf ace.passivation 
layer and the bulk even if no gate voltage is applied# and V is the applied gate 
voltage which we take with a coefficient less than unity# here for example with a 
weight of O.i. 

The calculation of surface recombination quantum i/f noise is similar to the 
calculation of quantum i/f noise from recombination in the space charge region. 
However, in this case the cross sections are not distributed over the width of the 
junction# but rather are concentrated at the surface which is caractenzed by the 
surface potential^. Therefore we can write an expression of the form 

Sjs(f) H e J s tanh x /f<T ho +7 po >exp(e$/kT>. (3.47) 

We note that the fabrication process of MIS structure introduces less bulk 
defects than the fabrication of photovoltaic devices. Nevertheless# the 
fabrication of MIS devices introduces some defects in the bulk layer located right 
under the surface. These defects will manifest themselves through a contribution 
to indirect tunneling. 
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4. Quantum i/-f Noise in -the Tunneling Rate 

In the case of tunneling -from a surface accumulation layer to the bulk, the 
velocity change of the carriers will lead us from the thermal velocity of the 
carrier on one side of the barrier to the thermal velocity of a carrier on the. other 
side of the barrier, if band to band tunneling is considered. If, however, the 
tunneling goes via intermediate states located in the bandgap, the velocity is 
zero as long as the carrier is stationary in the intermediate state. We can 
therefore write the i/f noise coefficient 

= (4ex/3'ip6kT/m*c^ (3.43) 

for band to band tunneling, and 

* (4*/3ffl3K7/m*c 2 (3.4?) 

for tunneling via intermediate states in the bandgap, where we have considered 
the average squared velocity change two times smaller. The effective mass is the 
mass of the minority carriers in the bulk material. 

For band to band tunneling from a surface inversion layer to the bulk, the 
velocity change of carriers corresponds to an energy difference of the order of 
the bandgap Eg plus an energy differnce of the order of the thermal energy 3kT/2, 
provided we are dealing with deep inversion, as used in practical MIS devices in 
the pulsed mode of operation. This yields the quantum i/f coefficient 

(4«/37T)(E g + 3KT/2)/m*c 2 (3.50) 

for band to band tunneling. For tunneling via intermediate states in the bancgap 
the corresponding energy differnce will be smaller, and therefore we replace E Q 
by Eg/2: 
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0^= (4*/3^r)(Eg/2 + 3KT/2)m*c 2 . (3.51) 

This relation is applicable both if the intermediate states are located in the 
depletion region or at the surface. 

In the last four equations we did not divide by the number of carriers 
simultaneously involved in the tunneling process, because this number is less 
than unity tor practical tunneling currents. Whenever a cross section or a 
process rate is tested with one electron or less than one electron at a time, the 
effective number of electrons in the denominator of the quantum 1/f formula must 
be replaced by unity. Let us calculate the average number of carriers 
simultaneously present in the tunneling process at any time. The tunneling 
process occurs over a distance d = Eg/eE, and the speed v of the carriers will be 
of the order of the thermal speed in the case of an accumulation layer, and of the 
order of the bandaap energy in the case of a surface inversion layer. Dividing d 
by v, we obtain a tunneling time t = 6'10“ 13 s for accumulation layers and t = 3- 
10" i3 s for inversion layers. From Eq. (2.8) we Know that the tunneling current is 
of the order of 10” 3 A/cm 2 . Multiplying this by t, we obtain 3 - 6*10“**’ C, i.e., 
2000 - 3000 electrons/cm 2 tunneling simultaneously in a device of 1 cm 2 area. In 
a device of dimensions 50 jm >: 50^.m = 2.5 10“ 5 cm 2 the average number of 
carriers in the process of tunneling at any time is therefore 0.05 - 0.075, and this 
is indeed much less than unity. Nevertheless, if the area of the device exceeds 5 
10“ 4 cm 2 , Eqs. (3.48) - (3.51) require an additional factor e/tJ*A which makes the 
noise spectral density proportional to and A, rather than to the square of 
these quantities. 

The photoelectric current will reflect the fluctuations in the number of 
photons arriving from the radiation background. The quantum efficiency will not 
exhibit considerable quantum 1/f noise, because the generated carriers will be 
corrected with certainty. Therefore the collection of photoelectrically generated 
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carriers is not controled by any cross section or process rate affected by 
considerable quantum 1 /f noise. If we neglect the i/f noise generated in the 
senes resistance of the diode, there should be no photoelectncally generated 1/f 
noise -from a short - circuited diode. 

Since the various dark current fluctuations with i/f spectrum are statistically 
independent, the total i/f noise is simply obtained by summing all contributions, 

IV. i/f NOISE LIMITED PERFORMANCE OF MIS DIODES 

From Eq. (2.4) we write the total dark current fluctuation in the form 

= ^Jdif + </J dep + '/ J s + / J tb + </ J tc + /' J tsc + (4,1) 

and the spectral density of current fluctuations will be neglecting/C 5 ^), 

s Jd = s Jdif + s Jdep + s Js + s Jtb + s Jtc + s Jtsc* lA ‘ 2) 

Here we have lumped the recombination current on the back surface J b together 
with the surface recombination (generation) current J £ . If we denote all the 
corresponding spectral densities of fractional fluctuations by a prime, S'ji = 
Sji/Ji 2 , we obtain 

S'jd = (t - T dif / ’V^dif + ^dep^d^Jdep + (J s /J d ,tS Js 
+ ‘Jtb 7 J d )2s 'jtb + < J tc /J d )2s 'jtc + <J t sc^d )2s 'jtsc* (4 - 3) 

This equation was obtained by dividing the previous equation through and 
shows that the biggest contribution will not necessarily come from the process 
with the highest fractional quantum i/f noise, i. e., with the highest i/f noise 
coefficient. The weight of each type of noise is determined by the corresponding 
squared current ratio. 
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Tl-e detectivity ot infrared detectors is limited in general by three types ct 
noise: u) curent noise in the detector, (Id noise due to background photons iphcton 
noise;, urn noise m the electronic system following the detector. We shall 
neglect here the background photon noise and the noise in the electronic system. 
The detectivity is defined as 

D*(X, t) = i,AA+; i/2 /MEP(cmHz 1/2 /W; i4.4) 

where A is the area Of the detector, NEP the noise equivalent power denned as 
the r.m.s. optical signal o-f wavelenghA required to produce r.m.s. noise voltage 
icurrent; equal to the r.m.s. noise voltage (current) in a bandwidth A 4, and 4 is the 
frequency of modulation. The noise equivalent power NEP is given by 

NEP = (hV/^q)( S Id <-fU*) 1/2 . ( 4 . 5 ) 

Therefore we obtain -for the detectivity 

D*<A, 4) = tyqfc/hc> CA/S Id (f)J 1/2 = O^qa/hc) CSj d <f>]' 1/2 (4.6) 

We notice that B*(a» f) is proportional to A up to the peak wavelength Ac. For A >A C 
we have 1 ^ = 0 and thus D*(7v, 4) = 0. By substituting our result tor Sj d , we 
obtain the general expression o 4 the detectivity as a -function of various 
parameters o-f the MIS device. 

Let us now evaluate the magnitude o-f the various dark current noise 
contributions. With m h * = 0.55 m D , m n * = 0.02 m D ,r n = 10“ 6 5, E g = 0.1 eV, 3kT/2 
= 0.01 eV, N s+ - = 10^/Vcm 2 , we obtain -for a p - type device wi+Ji Wp ]_ * 1 ' - 

S 'jdi44 = *Hnd + ^Hnr )Ce/ «n J diff :iF(a) sct coh <e 1/2 /C«iq^ 1/2 N p :}F(a)/a 
= (4.6 10" 3 /4fN p ) 4-10 _10 C 1/2 /C10~ 6 s. 1.5* 10 5 (cm 2 /Vs) 4-10 -21 J] 1/2 
= l.t'lO-^cm^/f, Cor = 10“* crn^/t, with 0( = 2 , 10 - ^ -for incoherent noise] (4.7) 

Up 

S 'jdep = C ^He e ^ (Zr no + ' C po ,J dep 3t4nh eV/2kT 
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= &C He e/feAwmeV/2kT)}eV/2kT = 0 *He /fAwn i = 4.6 10 _9 cm 2 /t (4.S) 

3 Js = <4<*/37r>(2/m*c 2 )C3kT/2 + eU/2 + 0.1 VeD He tanhx /^no + r po )J s 3 
■ (4^/37T0.02)(2/500000)C0.025 + 0.5 + 0.53 Ce tanhx /■feAwn i te eV/2kT - 1)] 

= 7 10~ 8 cm. 2 /f (4.9) 

= (4«#/37T)(E 0 + 3kT/2)/m*c 2 = <4/9.5 137 0.02H0.11 /500000) 

= 3.3 10“®cm 2 /t (4.10) 

S 'jlc = <^/3JT>(Eg + 3KT)/2m*c 2 « (4/9.5 137 0.02X0.12/10 6 ) 

= 1.8 10 _i:l cm 2 /f = S'j-tsc (4.11) 

Jdi-f-f was calculated in the small baas limit -for w p >>L, but w p = 0.25 L gives the 
same result; the incoherent case with a lattice constant of 0.65 nm and $= 320 K 
was also listed above (because a iOyun thick device is very short, so it may be 
applicable), and would give 1.8 10 -10 cm 2 /t tor a n - type device. Eqs. (4.8) - 
(4.11) would be reduced only m*/m R = 27.5 times tor n - type devices. We mention 
that S'j 5 has been calculated with the inclusion ot a term ot 10% ot the applied 
gate voltage V into the kinetic energy ot the carriers at the surface, and that tor 
the back surtace recombinaiton current this term has to be dropped in the similar 
expression ot S 'jb- However, we have neglected this here, because the surtace 
recombination terms will not turn out to be important, as we will see below. The 
applied gate voltage was taken to be V = 5 V. Using Eq. (4.3) and the current 
densities evaluated in Eqs. (4.10) - (4.14) to calculate the traction ot each 
current, we obtain 

1 cm -2 t S'j(t) = (20/132) 2 1.8 10" 6 + (10/132) 2 4.6 10~ 9 + (3.6/132) 2 7 10“ 8 + 
(0.01/132) 2 3.3 10“ 6 + (SO/132) 2 1.8 10~ 8 + (17.5/132) 2 1.8 10 -8 
= 3.67 10“ f; + 2.6 10" 11 + 5.2 10 -il * 1.9 10 -16 + 6.61 10“ 9 + 3.17 10" 10 
= 4.37 10 _s , or tor incoherent 1/t noise, 7.1 10 -9 (p) and 3 10“*® (n). (4.12) 
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This value can be used in order to estimate the detectivity of the device in our 
example. Substituting into Eq. (4.6), we obtain with a quantum efficiency 0.7 
and wavelength of iO^m 

D*a, f) = <^q^/hc) CS Jd (f)3 _1/2 = C0.7*1.6'1<T 19 C i0~ 5 m/(6.6'10“ 34 Js 3* 
10*m/s)j Cf/(4.37-10 -8 cm 2 i.74 > i0“ 6 A 2 /cm 4 )3 1/2 = 2*10 7 (cm Hz 1/2 /w> f 1/2 , or 
for incoherent 1/f noise, 5'10 7 (p>, and 2.5 10 c! <n>. (4.13) 

In conclusion we note that for the relatively large devices which we have 
considered, most of the quantum 1/f noise comes from fluctuations m diffusion 
and in the tunneling rate via impurity centers in the bandgap. The effective mass 
of the carriers is present in the denominator of all quantum i/f noise 
contributions except the coherent quantum i/f fluctuation present in the 
diffusion current of large devices. In smaller devices the diffusion current will 
also be given by the conventional quantum 1/f formula which contains the 
effective mass of the carriers in the denominator. For umklapp scattering the 
mass of the carriers in the denominator is even souared. Consequently we expect 
lower quantum 1/f noise from n - type devices, in which the minority carriers are 
holes, particularly if the devices are very small, e.g., below 10/tm. 

V. DISCUSSION 

The transition from coherent state quantum i/f noise to conventional quantum 
1/f noise is particularly interesting, and should be studied experimentally. This 
is possible with a sequence of devices Df smaller and smaller size, and will show 
a considerable change in noise at a size of the order of lowin'-. The theory of 
the transition is not yet well developed. Therefore, this experiment has 
particular importance: we do not know if the parameter s is sufficient to 
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characterize -the transission, and x-f -the parameter s should not be replaced by a 
power ofs, or by any other -function o-f s. The interpolation -formula used here is 
>jst a guess, or a speculation guided by the physical understanding o-f coherent 
quantum 1/f noise as a collective-field effect, and of conventional quantum 1/f 
noise as an effect which is not based on the collective field state of the 
particles, but arises from the individual field of each carrier. 

The most interesting component of the recombination current is the surface 
recombination current which plays a major role in the case of infrared detectors 
with pn junctions. Tn the case of MIS devices this role is not so important, as our 
calculation shows. Nevertheless, one should try to reduce both the concentration 
of recombination centers and the value of the surface potential jump U. This can 
be accomplished with careful surface treatment, and with a good passivation 
layer. Si02 layers have been successfully used by Radford and Jones in 
ion-implanted and double - layer epitaxial HgCdTe photodiodes^. 

In general the larger life time of the carriers in MIS devices, compared to 
junction devices is due to the absence of the damage inflicted by ion - 
implantation, or by the heavy doping required in double - layer epitaxial 
photodiodes. The quantum i/f noise is invesely proportional to this life time. 
Therefore, MIS devices should have lower 1/f noise. On the other hand, 1/f noise 
present in the applied gate voltage, in the timing of the readout and the value of 
the readout potential will be added as a 1/f noise source, if it is present. In the 
present calculation, however, this noise source has not been included. 

Any reouciion in the concentration of tunneling centers present in the bandgap 
will have a positive effect on quantum i/f noise. As we have seen in Sec. II, p - 
type devices should yield less tunneling via bandgap centers. The effective mass 
present in the denominator of the quantum i/f noise formula in this case should 
just be the effective mass of the carriers after the tunneling process, i.e., the 
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: fective mass of the outgoing carriers emerging -from the process we have 


considered, or the effective mass o-f the carriers coming in to the process of 


tunneling toward the centers in the bandgap. Here we have considered the 


tunneling process as the slower process which actually controls the rate of 


tunneling via bandgap states. The capture of carriers by the bandgap states is 


the second part of this compound process and has been considered fast enough, sc 


that it does not limit the rate of the total process. In general, however, both 


parts of the process have to be considered as a limitation on the rate, and in this 


case our noise formulae have to be revised through the inclusion of an additional 


term similar to tne recombination noise term. 


In the case of very small MIS devices, where only conventional quantum i/f 


noise should be present, we may find lower noise in the n - type devices, whose 


bulk minority carriers are holes with much larger effective masses than the 


electrons. This may happen in spite of the larger tunneling via bandgaD centers 


located right under the surface of these devices. 


Finally, we would like to empahsise that the present study has attempted to 


explain the basic concepts of quantum i/f noise and to illustrate their application 


to MIS infrared detectors. Although we have tried to pursue the calculation all 


the way to the evaluation of the detectivity, the data which we used in the 


calculation may not be applicable in the practical case at hand, and may have to be 


replaced with pertinent data in any concrete case. 


The author would like to acknowledge the help of M. Belasco, M. Kinch, E. Kelsi 


and R. Balcerak in many discussions on MIS devices and their noise problems. 









VII. APPENDIX B 


BACKGROUND OF THE G - R i/F NOISE DERIVATION 


At -first sight it looks as if (i) is incorrect since (R - G) should be proportional to 
the trap density Nj. As will be shown in the -following, the trap density is 
incorporated into the time constantsTp 0 andT n0 . 

According to Warner and Grung Cl4, eq. (4 - 191)3 

R - G « v p A p v n A n N t (pn - n^/v^^n + n £ ) + v p A p (p + n*) (B 1) 


Here, v p is the average speed o-f the holes, v p the average speed -for electrons, 
A p the trap cross section -for holes, and A n the trap cross section -for electrons. 


Multiplying denominator and numerator in (3 1) by N^ and then dividing both by 
v o A o N t v n A n N ty ield5 
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R - G = (pn - n^j)/C<n + r^rr^ + (p + 
which is just Eq. (3.23). in addition 

T nn = i/v p A p N t = i/C D , C n = v D A D N t 


po 


i/v n A n N t = 1/C n , 


C n = v n A n N t 


as had to be proved. We see thatTr po and ^ nD are real - time constants 


1 : c;” 1 ). 


process. According to Van Vliet £153 the true g - r time constant "T^ is 

= £<n + nj)v n A n + <p + n^VpAp] = N^/C(n + n^C,^ + (p + n^CpD (B 3) 

It we put i/T n = <n + nj)C n /N t , l/r p = (p + n^Cp/N^. then 

T t* 7 pV'V r n>- s3 *> 

Our Eq. (3.23) is thus tully correct. It N* goes to zero,T no and Tp 0 go to 
intinity and hence (R - G) goes to zero. Note, however# that 2"^ is independent ot 
N^i except tor. a possible dependence ot n and p upon N^. 

'We also note that A p and A n are practically equal. For the case that the trap 
is negatively charged when occupied by an electron, the excited electron interacts 
with a neutral trap and hence has a cross section ot the order ot l(T l6 cm“. But 
when a hole is created by having an electron trom the valence band recombine with 
the empty trap, the electron also interacts with a neutral trap and hence has 
practically the same cross section. It the trap is neutral while occupied by an 
electron, the excited electron interacts with a positively charged trap and hence 
has a cross section ot the order ot 10 - ^ - 10“^cm‘-. But the same is true when 








an electron from the valence band interacts with a positively charged trap. 
Hence, and A n are always nearly equal. 

However, Vp is inversely proportional to unp*)*^ 2 anc j Vn is inversely 
proportional to <m n *) i/2 , so thatT no is proportional to <m n *> 1/2 and to 
(m_*)^ 2 . Consequently 


V^^^no 1 ^ m */C(m n *> 1/ - * <m 


We need this for the calculation of°^ ^ in the g - r i/f noise theory. 

There are two forms of quantum 1/f noise. In the first place Cp and C n 
fluctuate because A p and A n fluctuate independently in a i/f fashion, according 
to quantum 1/f noise theory. This would yield Hooge parameters eC' H and 


dt' H p = <A<K/3 71 )Uv p 2 /c 2 ) = (4^/37n(3kT/m p V) 




m 


et' Hn = <4^/37t)(4v n 2 /c 2 > = (4^/3Jn(3KT/m n # c 2 ), 


since rrip fr AVp 2 /2 = m n *Av n 2 /2 = 3/2kT for the excitation of electrons and holes 


In the second place the electrons and holes in the space charge region are 
decelerated or accelerated. The potential across the space charge region of width 
w'p.f.f is a(V d ^ - V), where a = 1/2 for a symmetric Unction and aC^ 0.75 fo an n + - 
p junction. Electrons excited from the trap level are accelerated and thus gain an 
energy E n = ea(V cj j^ - V) = m n *Av n 2 /2 and holes excited from the trap level are 
accelerated and gain an energy E p * e(l - - V) = mpVvp^/Z. This 

produces quantum 1/f noise such that 










-so 


^Hp = C2e(l - a)(V dif - V)/m p *c 2 3(4ol/37l) (B 5) 

cl"Hn = C2ea(V dif - V)/m n *c 2 3(4«</37r> (B 5a) 

Adding -the two independent contributions toot Hn and^ Hp yields 

Hn = <4«</37r){C2ea<V dif - V + 3kT3/m n *c 2 > 

d Hp = <4^/37TKC2e(l - a)(V dif - V) + 3kT3/m n *c 2 }, (B 6) 

so that -from (B 4), (B 6), and (3.35) 

H * Wpo'ttpo t ~no |2 *Hp * £ 7 no' ir po +Tno» 2 *Hn 
■ (4oi/3nKC2e(V dif - V) + 6kT3/C(m n *) 1/2 + <m p *) 1/2 3 2 c 2 . (B 6a) 

For Hg^j.CdjjTe with x = 0.20, m n * = O.Oim and m p * = 0.55m so that C(m n *)^ 2 + 
(m p # )* /2 3 2 = 0.7im which is a factor 7i larger than m* n . An earlier estimate by 
Radford and Jones C133 used m n * instead of C(m n *) 1/2 + <m p *) 1/,2 3 2 . 
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